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Abstract 



This paper aims to review our recent results on exploring the capabilities of nonquantum field 
theory as a possible tool for describing single photon-like objects, considered as massless time- 
stable spatially finite physical entities with compatible translational-rotational dynamical struc- 
ture. It consists of five chapters and includes 22 sections and subsections. 

In Chapter 1 (Introduction) we present briefly some remarks on Maxwell vacuum equations 
and our notions concerning the two basic concepts in physics: physical object and interaction, 
with an accent on the view that physical interaction necessarily implies energy-momentum ex- 
change, and on the suggestion that the energy-momentum exchanging subsystems of a general 
field should be described rather by (F, *F), than by (E, B). 

Chapter 2 (Nonrelativistic considerations) begins with a reasoning on the status, right 
understanding and appropriate use of the Coulomb force law in electrostatics, leading to the 
conclusion that the usual way of introducing static electric field as local object directly from 
the Coulomb force law violates in definite sense the local conservation laws in the frame of 
Maxwell electrodynamics. Then we concentrate on the role, significance and eigen properties 
of the Maxwell stress tensor in nonrelativistic terms. This part of the paper culminates in 
writing down nonlinear field equations for the vacuum electromagnetic fields paying due respect 
of the Newton view on the sense of dynamical equations as local balance relations of conserved 
quantities carried by two appropriately defined subsystems of the general electromagnetic field, 
and presenting some important properties of the nonlinear solutions. Finally, we discuss briefly 
the contents of this part of the paper. 

Chapter 3 (Relativistic considerations) of the paper makes use of the relativistic for- 
malism. We begin with presenting the notion of photon-like object(s) (PhLO). The existing 
corresponding relativistic description of PhLO in the frame of Extended Electrodynamics is 
briefly recalled. In terms of integrability and nonintegrability properties of distributions (dif- 
ferential systems) on a manifold a principle of geometric interaction between two nonintegrable 
distributions is formulated and a corresponding physical interpretation is given. This principle 
is further substantially used in building mathematical description of single PhLO, considered 
as composed of two individualized and interacting subsystems, and realizing a special kind of 
dynamical equilibrium. The mathematical model is built on the assumptions that Frobcnius 
integrability can be made to correspond to physical time-stability, and the nonintegrability of 
sub distributions of integrable distribution to correspond to local physical interaction between 
subsystems of a time-stable continuous physical system. Two approaches are considered: direct 
use of the Frobcnius theorem and the corresponding curvature being a measure of nonintegra- 
bility, and the recently developed approach known as "non-linear connections". Both these 
approaches make use of the corresponding curvatures for generating appropriate quantities de- 
scribing local physical interaction, i.e. local energy-momentum exchange 

A corresponding concept of electromagnetic strain is defined and the basic stress-energy- 
momentum relations, obtained before, are represented in electromagnetic strain terms. 

Chapter 4 (Equations of motion for PhLO. Solutions) deals with various views on dy- 
namical equations for PhLO and gives illustrative examples of appropriate solutions. 

Finally, Chapter 5 (Retrospect) gives a retrospective discussion of the contents of the paper. 
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Chapter 1 

Introduction 



1.1 Preliminary remarks 

As it is well known the vacuum Maxwell equations (zero charge density: p = 0) do not ad- 
mit spatially finite time-stable solutions of photon-like type, i.e. solutions, having compatible 
translational-rotational dynamical structure and propagating as a whole along straight lines in 
the space with the fundamental velocity "c" - the speed of light in vacuum and without dis- 
persion. This is in corressponding degree due to the fact that in the frame of Maxwell vacuum 
equations every component U{x, y, z, t) of the electric E and magnetic B fields, satisfying cor- 
responding smoothness conditions, necessarily satisfies the D'Alembert wave equation UU = 0, 
and according to the Poisson's theorem for this equation, every spatially finite initial condition 
J7(a;, y, z, 0) = (p(x,y,z); ^{x,y, z,0) = ^{x,y,z), where (p and ^ are finite functions satisfy- 
ing definite differentiability conditions, blows up radially and goes to infinity with the speed of 
light [1,2,3]. So, in such a case, through every spatial point outside such an initial condition 
pass fore-front and back-front, and after this the corresponding point forgets about what has 
happened. Hence, photon-like objects require new kind of description. 

On the other hand the more than a century successful usage of Maxwell equations in de- 
scribing various electromagnetic physical systems and processes undoubtedly suggests that their 
adequateness to reality is hardly understood and used fully. For example, one of the crucially 
important relations that follows from Maxwell equations, and that has been successfully used 
throughout all these years, is the Poynting energy-momentum balance equation in vacuum 



where c is the velocity of light in vacuum. Considered from the above mentioned point of view it 
turns out that this Poynting equation admits time-stable, spatially finite and propagating along 
straight lines solutions of the following kind 

E = [u{x, y, ct + £z), p {x, y, ct + ez),0], B = [ep (x, y, ct + ez), —£u{x, y, ct + ez),0], e = ±1 , 

where u and p are arbitrary functions, so they can be chosen finite. This observation suggests 
to look deeper and more carefully into the structures and assumptions used for mathematical 
interpretation of the experimental electric-magnetic induction discoveries made in the 19th cen- 
tury. In other words, which relations, and on what grounds, should be defined as basic, so that 
the further deduced equations and relations to give reasonable and physically acceptable results 
when viewed from the viewpoint for spatial finiteness and time stability. Finding the right way 
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to choosing adequate mathematical objects and corresponding equations seems specially im- 
portant when we try to describe the intrinsic dynamical properties of such spatially finite and 
time stable field objects. Therefore, it seems very important to have the right notion about 
concepts like physical object, intrinsic property, dynamical property, identifying characteristics, 
dynamical structure, admissible changes, field equations, etc. 

The idea to extend the vacuum Maxwell equations in such a way that spatially finite time 
stable and straight-line propagating solutions to be incorporated is, of course, not new [4]. 
Moreover, a general principle concerning all theoretical physics was formulated by Born and 
Infeld [5] , stating: a satisfactory theory should avoid letting physical quantities become infinite. 
Other recent efforts in this direction in the frame of electrodynamics may be found in [6,7,8,9]. 
We are not going to analyze here all these various approaches, what however deserves to be 
noted is that the new equations offered therein have not usually direct physical sense of energy- 
momentum balance relations as the Newton law in mechanics has. And this is true also for 
the very Maxwell equations. The above mentioned example with the Poynting relation clearly 
suggests a more serious and physically motivated respect to be paid to this aspect of the problem 
when some basic relation in a physical theory is postulated. For example, should we consider 
the very electric and magnetic fields (E,B) as force fields in case of no charges present, or 
the divergence VjM*-' of Maxwell's stress tensor M*-' [10] should be accepted in the pure field 
case as force field? Each of these three vector fields generates integral lines, so, which integral 
lines should be considered as force-lines, and what is a force-line in case of absence of charged 
particles? Clearly, since the general concept of force, considered as local object, is local energy- 
momentum exchange, more reliable seems to be the divergence ViM^^, so, the Newton-likc 
equations of motion in the vacuum case should read V^M*-' = 0. In other words, the implied 
energy-momentum exchange between the electric and magnetic components of the field should 
be represented by relations having direct energy-momentum exchange sense, and not by the 
local versions of the so called induction laws 

dB ^ ^ dE 

rotE+ — = 0, rotB- — = 0, 

where ^ = ct. In fact, we can NOT observe and verify directly these relations since we have 
not corresponding devices, we verify them indirectly through observing corresponding energy- 
momentum changes, i.e. state-behavior changes, of charged particles. In other words, we have a 
system having two subsystems: field and charged particles, and, assuming the energy-momentum 
local conservation law, wc make energy-momentum change conclusions about the field through 
the corresponding mechanical energy-momentum changes of the charged particles. Therefore, 
when we want to understand and describe intra-field energy-momentum exchanges, i.e. in case 
of no charged objects presented, we must have in mind some preliminary pictures about the 
structure and possible more or less individualized and time-stable substructures among which 
this energy-momentum inter-exchange should take place. Assuming the electric and magnetic 
components as such substructures we find that the absence of well defined local interaction energy 
between these two components (recall that the energy density of the field is ^(E^-)-B^)) seriously 
complicates this problem: the two Maxwell equations above imply some energy exchange between 
E and B, so how it is performed if there is NO local interaction energy, moreover, neither E 
nor B are capable to carry separately momentum and angular momentum, so, which are the 
corresponding energy-momentum exchanging substructures of the field during propagation? 

The relativistic development of classical electrodynamics introduced new point of view: ad- 
equate mathematical objects that represent such two substructures of the general vacuum field 
are not E and B, but two differential 2-forms -F'(e,b) and *-F(_b,e) on Minkowski space-time, so, 
from the new point of view, any internal energy-momentum exchange should take place between 
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F and *F. However, the "new" field equations dF = 0, d * F = 0, although in terms of F 
and *F, keep the old viewpoint, and in a definite sense they forbid such internal exchange (see 
Sec.3.1). This study is directed to find corresponding new equations. 

The paper is organized as follows. In Sec.1.2 we present some general considerations concern- 
ing the mathematical description of a physical system. Sections 2.1-2.5 are devoted to defining 
the model in pre-relativistic terms. Sections 3.1-3.6 give relativistic approach, making use of 
Frobenius curvature, non- linear connections and relativistic strain. Sections 4.1-4.4 present var- 
ious views on dynamical equations and give solutions of photon-like nature. Finally, Ch. 5 gives 
a retrospect of this study. 

1.2 Physical Objects and Interactions 

From a definite point of view every physical system is characterized by two kinds of properties. 
The first kind of properties we call identifying, they identify the system throughout its exis- 
tence in time, so, the corresponding physical quantities/relations must show definite conserva- 
tion/constancy properties (with respect to the identification procedure assumed). Without such 
experimentally established properties we could not talk about physical objects/systems at all. 
The second kind of properties (which may be called kinematicat) characterize the time-evolution 
of the system, the corresponding quantities are time-dependent, and the corresponding evolution 
is consistent with the conservative/constant character of the identifying properties/quantities. 
In this sense, the equations of motion of a physical system can be considered as relations deter- 
mining the admissible time-changes of these quantities. For example, the mass m of a classical 
particle is an identifying quantity, while the velocity v is a kinematical one. This view implies, 
of course, that the external world acts on the system under consideration also in an admissible 
way, i.e. an assumption is made that the system survives, the interaction with the outside world 
does not lead to its destruction. 

In theoretical physics we usually make use of quantities which arc functions of the identifying 
and of the kinematical characteristics of the system and call them dynamical quantities. A 
well known example is the momentum p of a particle: p = mv. Of crucial importance for the 
description of admissible changes are the conservative dynamical quantities, i.e. those which may 
pass from one physical system to another with NO LOSS. For example energy and momentum 
are such quantities, moreover, they are universal in the sense that every physical object carries 
non-zero energy-momentum and, vice versa, every quantity of energy-momentum is carried by 
some physical object. So, if a definite quantity of energy-momentum passes from one object to 
another, this same quantity of energy- momentum can be expressed in terms of the characteristics 
of the two objects/systems, and the two expressions to be equalized. This allows to describe 
interaction between, or among, physical objects. Thus we have a consistent with the requirement 
for "identification through conservation" -way to write equations of motion, and this is the way 
used by Newton to write down his famous equations p = F, where F carries information about 
where the momentum change of the particle has gone, or has come from. This also clarifies 
the physical sense of the concept of force as a change of momentum, or as a change of energy- 
momentum in relativistic terms. Paying due respect to Newton we shall call some equations of 
motion of Newton type if on the two sides of " =" stay physical quantities of energy-momentum 
change, or energy-momentum density change in the case of continuous systems. Note that, 
written down for the vector field p, i.e. in terms of partial derivatives, the above Newton 
equation looks like VpP = mF, where the left hand side means performing two steps: first, 
determining the "change quantity" Vp, second, projecting Vp on p, and the right hand side 
may be expressed as a function of the characteristics of both: the particle and the external 
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physical environment. 

If there is no energy-momentum (or energy-momentum density) change, then putting the 
corresponding expression equal to zero, e.g. VpP = 0, we obtain the "free particle" or "free 
field" equations. In such a case we just declare that only those changes are admissible which are 
consistent with the (local and integral) energy-momentum conservation. 

We note that an initial extent of knowledge about the system we are going to describe 
mathematically is presupposed to be available, so that the assumptions made to be, more or 
less, well grounded. This knowledge is the base that generates corresponding insight and directs 
our attention to the appropriate mathematical structures. This is exclusively important when 
we deal with continuous, or field, physical objects/systems. 

To illustrate our idea, let's consider a many-component continuous system, i.e. such that 
each component is assumed to be time-stable and recognizable during the system's existence. 
Accordingly, the wholeness-|-structural integrity of the system should be mathematically 
represented by some complex <I> of interconnected fields: $ = {^a},CL = 1,2,..., where each 

represents some individualized subsystem and may be also many-component one: = 
. . .}. In view of the above considerations, if our system is free, the steps to follow are: 

1. Specify and consider the mathematical model-object $ = . . . }, a = 1, 2, . . . which 
is chosen to represent the wholeness-|-structural integrity of the physical system considered; 

2. Define the change-objects !?($„), considered as due to internal interactions; 

3. "Project" D{^a) on ^a/^b:^ ^1 by means of some (in most cases bilinear) map *P; 

4. The projections ^{D(^a),^b) and *P(£)($5), $„) obtained, we interpret physically as 
local energy-momentum exchange between the individualized subsystems described by and 
^b- V{D{^a),^b) is the energy-momentum that transfers to and ^{D{^^),^a) is the 
energy-momentum that $5 transfers to 

5. The subsystem described by iiiay, or may not, keep its energy-momentum unchanged 
during the internal interaction. Correspondingly, in the first case we'll have ^{D{^a), ^a) = 0, 
and in the second case we shall have ^{D{^a)i^a) 7^ 0- 

6. The sum Sa^P(I?($a), <I'a) should be equal to zero, meaning that our system $ = {$a} 
conserves locally its energy-momentum. 

The zero value of the projection ^{D{^a), *^a) is interpreted in the sense that the identifying 
characteristics of have not been disturbed, or, the change -D($a) is qualified as admissible. 
This consideration shows the importance of knowing how much and in what way(s) a given 
physical system is potentially able to lose, or gain energy- momentum (locally or globally) , without 
losing its identity. 

It is always very important to take care of the physical sense of the quantities that we put on 
the two sides of the relation A = B. Mathematically, from set theory point of view [13], A and 
B denote the same element, which element may be expressed in different terms, e.g. the real 
number 2 can be expressed as 3 — 1 = 6/3 = 7^ (2.x + const) and also in many other ways. From 
physical point of view, however, we must be preliminary sure that A and B denote the same 
thing qualitatively and quantitatively, i.e. A and B must denote the same physical quantity. This 
is specially important when the equation we want to write down constitutes some basic relation. 
And the point is not the physical dimension of the two sides to be the same: any two quantities 
by means of an appropriate constant can be made of the same physical dimension, but this is a 
formal step. The point is that the physical nature of the physical quantity on the two sides must 
be the same, and this should be well understood and correspondingly guaranteed beforehand. 

For example, it is quite clear that on the two sides of the Newton's law p = F stays the well 
defined for any physical system quantity " change of momentum" since the momentum quantity 
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is a universal one. For a counterexample, which physical quantity stays on the two sides of the 
Poisson equation A?7 = kp, k = const? On one hand, such a quantity is expressed through AU 
and, since grad U is usually interpreted as force, AU appears as a "change of force" characteristic 
of the field U since it is essentially defined by the second derivatives of U. On the other hand, the 
same quantity is expressed through kp and appears as a characteristic of the mass particles, so, 
do we know such a quantity? The same question can be raised for one of the Maxwell equations: 
rotB- iE = ^j. 

c c •' 

In the case of classical particles momentum is always represented as the product mv and this 
is carried to fluid mechanics (continuous mass distribution) as p,{x, y, z; t).v{x, y, z; t), where jj, is 
the invariant mass density. A similar quantity is introduced in electrodynamics as electric current 
density j = p{x, y, z; t).-v{x, y, z; t), where p is the electric charge density. The energy-momentum 
exchange between the field energy-momentum and the mechanical energy-momentum of the 
available charged particles is described by the force field F = pE-l- ij x B. So, the corresponding 
Faraday-Maxwell force lines should be the integral lines of the vector field F. Clearly, in the 
charge-free case we get F = 0, so the concept of force-lines defined by F does not work. Hence, 
if we would like to use this concept appropriately in the charge-free case, we have to introduce 
it appropriately. The simplest way seems to consider the integral lines of E and B as force lines 
also in the charge free case, but we do not share this view: if p = then j = 0, the force-vector 
is zero and NO integral force lines exist. The two vectors E and B generate, of course, integral 
lines, but these integral lines are NOT force lines in the vacuum case since the vacuum E and 
B are NOT force fields, and such an interpretation of the integral lines of E and B would be 
misleading. In fact, if in case of E some compromise could be made since E and pE are colincar, 
in the case of B this is impossible: the corresponding force lines are generated by j x B, so, at 
each point they are orthogonal to the integral lines of B. We note that these problems arise 
in connection with passing from discrete (point-like) quantities to continuous quantities: mass 
and charge distributions, currents, etc., so such a transition must be carefully performed in any 
theory in view of the above remarks. The rule we are going to follow is: the sense of quantities 
used must always be quite clear and must not be misleading. 

As we already mentioned, there exists a sufficiently good force field defined by Maxwell in 
terms of the divergence of his stress tensor M*-^ , which definition works quite well also out 
of and away from any media built of, or containing, charged mass particles. So, 
in the frame of the theory at the end of 19th century if we ask the question: if there are NO 
charged particles and the time- dependent EM- field cannot transfer energy-momentum to them 
by means of the force field F, and the propagation of the free EM-field is available, so that 
energy-momentum internal exchanges should necessarily take place, how these processes and the 
entire propagational behaviour of the field could be understood and modeled^., the right answer 
in our view should be: turn to M*^ and consider carefully the divergence VjM'^ terms 
as possible force fields generating corresponding force lines along which energy- 
momentum is internally and locally transported between/among subsystems. As 
will be seen further in the paper, such a look on the issue would necessarily lead Maxwell 
and his followers to the prediction that real, free, spatially finite and time-stable formations 
of electromagnetic field nature having compatible translational-rotational dynamical structure 
should exist, a result that has been proved in studying the photoeffect phenomena about 30 
years after Maxwell's death. 

We consider as a remarkable achievement of Maxwell the determination of the correct ex- 
pressions for the energy density of the electromagnetic field through the concept of stress [10]. 
His electromagnetic stress tensor IVP^ still plays an essential role in modern electromagnetic 
theory as a part of the modern relativistic stress-energy-momentum tensor. However, by some 
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reasons, Maxwell did not make further use of the computed by him divergence ViM*-' of the 
stress tensor (and called by him "force field" [10]) for writing down Newton type equations 
of motion for a free electromagnetic field through equalizing different expressions for the same 
momentum change. Probably, he had missed an appropriate interpretation of the vector c E x B 
(introduced by Poynting 5 years after his death and called "electromagnetic energy flux" [14]). 

In connection with the above considerations the following more general question arises: docs 
theoretical physics make the right step allowing static force fields to be written on the right hand 
side of Newton equation of motion in mechanics? Every static field, whatever is its nature and 
origin, necessarily conserves locally and globally all its physically meaningful characteristics, so 
its energy is also conserved and is not at disposal to other physical systems. Moreover, static 
fields do NOT propagate, so they do NOT carry momentum, therefore, there is NO WAY other 
physical systems to lose or gain momentum at the expense of externel static fields. Hence, trying 
to pay respect to the momentum conservation through writing down Newton type dynamical 
equations p = F with static force field F, we, in fact, violate it: due respect requires due usage, 
so, in our view, no selfrespecting theory should allow dynamical equations with static F! 

This concerns not only mechanics, and in fact, not only static fields: propagating compos- 
ite time-stable physical systems are supposed to consist of individualized time-stable subsystems 
capable to interact, i.e. to carry and exchange energy-momentum during propagation. Do we 
respect always this natural principle in physical theories? Consider for example the two vector 
vacuum Maxwell equations: B = — crot E, E = crot B. These equations imply that there is mu- 
tual physical influence between E and B, which we understand physically as energy-momentum 
exchange between the individualized E and B components: E acts upon B and B acts upon 
E. Now, according to the equations, each of these two components propagates and keeps its 
individualization during propagation, so, it should be able to carry momentum. However, the 
energy concept of the theory excludes nonzero interaction energy between these two compo- 
nents to exist since the energy density w is given by the sum of the energies carried by E and 
B: -u; = ^(E^ + B^), so, how does energy exchange take place? Further, the momentum concept 
in the theory is defined and experimentally proved quantitatively to be given by ^(E x B), 
so neither of the assumed in the theory electric and magnetic components is allowed to carry 
momentum separately. In view of this, how e.g. in the plane wave solutions, where the re- 
lations E^ = B'^ and E^ + B^ = 2|E x B| always hold, the implied by the above equations 
internal energy-momentum exchange between the two subsystems mathematically identified as 
E and B is performed? May be we have not made the right mathematical identification of the 
subsystems? 

The suggestion that we should come to is that we must very carefully make significant conclu- 
sions and assumptions about the dynamical structure of the physical system under consideration, 
especially in case of continuous and spatially propagating composite systems. In this work we 
shall try to follow the rule that an adequate choice of mathematical structures in the theory must 
correspond to a sufftciantly well recognized, identified and interacting physical structures. In the 
spirit of this, our hope is that, in general, an isolated time-stable nonstatic continuous physical 
system should correspond to a completely integrable distribution A on an appropriately defined 
manifold, its interacting subsystems Ei, E2, ... should correspond to NONintegrable subdistribu- 
tions Ai, A2, ... of A, and the very interaction, i.e. local energy-momentum exchange, between 
any two subsystems Sj <^ T,j can be described mathematically in terms of the corresponding 
Frobenius curvatures ria. 
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Chapter 2 

Nonrelativistic Considerations 



2.1 How to understand the Coulomb force ? 

Usually, the Coulomb force field is introduced starting with the Coulomb force law: / = 
where q and Q are the charges of two small bodies (usually considered as point-particles) and 
r is the euclidean distance between them. Two fields Eq and E^, considered as generated 
correspondingly by any of the two charges Q and q, are defined by the relations 



Now the so defined quantities Eq and Eg are considered as vector fields, i.e. local objects, 
defined outside the regions occupied by the source objects of charge magnitudes Q and q, and 
are interpreted as force-fields acting on other unit charges, hence, the force acting on the ^'-charge 
is gEg and the force acting on the Q-charge is QEg. 

The mechanical behaviour of the q'— particle in the reference frame connected with the Q- 
particle is defined by the Newton law Pg = qEq, where Pg is the mechanical momentum of 
the g-particle in this reference frame (clearly, pQ = 0) . This " dynamical" equation presupposes 
that the change of the mechanical momentum of the g'-particle comes from (or goes to) the 
corresponding change of the momentum carried by the field Eg in accordance with the universal 
momentum conservation law. However, such a justification assumes that the field Eg carries 
non-zero energy and momentum and is capable to exchange them with other physical objects. 
How much is, for example, the field momentum? This question requires corresponding definition 
of the field momentum. This directs our attention to the theory based on Maxwell vacuum 
equations. Maxwell theory, however, gives three objections to this understanding of the physical 
situation: 

-the field EQ(a;, y, z), considered as local object, i.e. vector field, satisfies the vacuum equa- 
tions rot Eg = 0, div Eg = outside its source, and according to the theory, every such solution 
field conserves its energy, momentum and angular momentum; 

-the static nature of the field requires no time-changes of any field characteristic; 

-the field momentum density in the theory is proportional to the Poynting vector, so, neither 
the electric field E nor the magnetic field B are allowed to carry momentum separately. 

In general, every vacuum solution of Maxwell equations conserves its energy, momentum and 
angular momentum, so, NO vacuum solution (E, B) should be allowed to participate 
directly as force generating agent in the expression gE + |v x B. 

We see that for vacuum fields the usual setting "charged particle in external field" does not 
work: the field can not afford any chance to the "test particle" to gain locally momentum from 
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the field. In view of this, how to understand the Coulomb force law from theoretical point of 
view in the frame of Maxwell theory? 

In order to answer this question we make the following considerations. First, some clarifica- 
tions concerning the structure and admissible changes of the physical situation. We have two 
mass particles carrying electric charges q and Q. The two masses "generate" two gravitational 
fields which are further neglected as physical factors. The two charges "generate" two electric 
fields: Uq - denoted further just by u;, and fig - denoted further just by Q,. The whole system is 
static and time stable, so, the two fields and the two particles considered as mechanical objects, 
exist consistently with each other. 

Remark. We put the term "generate" in commas intentionally, because in this case the charge- 
field configuration we consider as the real one, i.e. no charged particle can exist without such a 
field, and no such a field can exist without charged particle. In the theory the opposite of the 
usually stated idea that charges generate fields is realized: not the charged particle "generates" 
field, but the field "generates" charge through the Gauss theorem, so, both the charge and the 
field aspects of the situation should be paid equal respect. 

Since we consider electrostatic situation, no magnetic fields are assumed to be present. The 
admissible changes, by assumption, do NOT lead to destruction of any of the objects. Paying 
now due respect to the Gauss theorem we have to assume that each of the two fields is NOT 
defined inside the small region that its source occupies. Therefore, the two fields u and Q, are 
defined on the topologically non-trivial space S = M?\{WqUWQ), where Wq and Wq are the two 
small nonintersecting regions, treated further as two balls with boundaries 5^ and Sq, occupied 
by the two particles. How to specify the mathematical nature of u and Q, ? 

The topology of S, which must be kept unchanged, and the assumed spherical symmetry 
of each of the two fields with respect to (the centers of) Wq and Wq suggest to choose the fields 
n and oj as spherically symmetric representatives of the 2-dimensional cohomology group of S. 
We introduce two spherical coordinate systems {r,9,(p) and {f,9,(p), originating at the centers 
of Wq and Wq respectively, so, any two spherically symmetric with respect to the centers of Wq 
and Wq 2-forms will look as follows: 

Q{r, 6, ip) = h{r) sin 9d6 A dip, uj{f, 9, ip) = h{r) sin 9d9 A dCp. 

Being representatives of corresponding cohomology classes and lv must satisfy dO. = and 
da; = 0, so, h{r) = Const and h{f) = const. We denote Const = Q and const = q. Now, the 
euclidean Hodge star operator * and the euclidean identification of vectors and covectors give 

Eq = *0(r, 9, ip) = ^dr, = *u{f, 6, (p) = ^df. 

Going further we note that at every point of E the real field configuration is built of two 
physical fields of the same physical nature, therefore, the resulted stress should depend on the 
local mutual interference/interaction between the two stress generating fields and u. The 
point is how to model mathematically this local interaction of the two fields? At this moment 
the Maxwell stress tensor M/ = Z^Z-^' - \Z'^6{ = \n{Z)rnn n{Z)"'''6{ - n{Z)im Vt{Zy'^, 0(Z) = 
*(Z), which is defined by any vector field Z, could help us as follows. 

Mathematically, the tensor M can be considered as a quadratic map from the vector fields 
on S to (1, l)-tensors, i.e. to the linear maps in the linear space of vector fields. Each of our two 
fields generates such (1, l)-tcnsor field: M{uj) and M(J7). Recall now that every quadratic map 
$ between two linear spaces generates a bilinear map r$ according to T$(x,y) = <I>(x + y) — 
— where (x, y) are corresponding variables. So, in our case we can define corresponding 
bilinear map. 
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Identifying the vector fields and 1-forms on S by means of the euclidean metric e: E' = e*^ Ej, 
for the two Maxwell stress tensors, expressing here stresses of topological origin, we have 

Mq = M(Eq) = Eq ® Eq - ^ E^ id^s, Mq = M(Eq) = Eq ® Eq - ^ E^ iciTS. 
The corresponding bihnear map T will be 

T(Eg, Eg) = Eg ® Eq + Eq (8) Eg - Eq-E^ z^ts- (2.1) 
In components we have correspondingly 

T^E^Eq) = {E.UEqY + {EQUEgY-Eg.EQSl 

The tensor field (— ;^T) may be called mutual stress tensor, or interaction stress tensor. In 
r*S(g)rS we have the trace form tr, and on S we have the standard volume form Ug = dxAdyAdz. 
So we can form the object tr (8) cOg- By definition, the quantities 



w 
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{tr(SiLOo){-—T) = -—<tr,T>u;o, and U= w (2.2) 
47r 47r 7s 



will be called interaction energy density and interaction energy for uj and ft. Clearly, w and 
U may be positive, zero, or negative. Further we shall follow the rule that an isolated (qua- 
sistatic) physical system of this kind tends to configurations with minimal value Umin of the 
integral interaction energy U, hence, an intrinsically induced drifting between two allowed static 
configurations should satisfy the relation SU < . 

In order to compute U we compute first w and obtain 



1-1 1 
w = — E.Ea;o = — A *a; + a; A *0) = - — 
Att Stt Stt 



Making use of the Stokes theorem the integral of w over E is transformed to 2-dimensional 
surface integral over the boundary dT, of S: 5S = S^^^^ U S^U Sq. On S*^ the corresponding 
integral has zero value. So, in the induced on dT, orientation, and denoting by Rq and Rq the 
radiuses of Sg and Sq respectively, we have 

^ q 1 f Rl^QsinOdO A d(p Q 1 f R^qsmBdO Adip 



r RqQ sm.U dO A d(p Q I r 
2 AttRI Jsivjsl r ^2 AttR^ Jg^^sl 



''Q 

On Sq we have f = const and J^2 = 0. Similarly, on Sq we have r = const and Jg2 oj = 0. 

Notice further that Ms a harmonic function, so, at every point p ^ T, it can be represented by 
its average value on the corresponding 2-sphere centered at p. Now, the first integral reduces 
to integral over the 2-sphere Sq and it is equal to similarly, the second integral reduces to 

integral over the 2-sphere 5^ and has the same value, where R is the euclidean distance 
between the centers of the two spheres. Thus, the computation gives finally U = 

Now, according to the above mentioned rule that 6U < 0, and that q and Q do not change, 
for the case q.Q > we obtain SU = —^SR < 0, so SR > 0, i.e. repulsion should be expected; 
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and for the case q.Q < we obtain SU = —^6R < 0, so 6R < 0, i.e. atraction should be 
expected. 

The above consideration clearly suggests the idea that the Coulomb force law originates 
from available interaction between the two fields Eq and Eg under quasistatic changes of the 
integral interaction energy U leading to minimization of U. In fact, if U changes then the change 
SU = —^5R must be carried away mechanically by the (g, m)-particle: = 5U, since there 
is no other physical factor in the system considered. So, the Coulomb force can be understood 
as an integral characteristic of the system, therefore its field, i.e local, interpretation may be 
reconsidered. On the other hand, in the corresponding spherical coordinates, {Q * a;) and {q * CI) 
look very much as SU, but this first-sight resemblance should not mislead us. The difference 
is quite serious: (Q * u) and {q * 0,) are 1-forms, local objects by definition, while ^SR is not 
local object, R is not the coordinate r and, contrary to dr, SR is not 1-form on S: we should 
not try to obtain local objects just through noting some possible change tendencies of integral 
characteristics of a system. We may allow ourselves to call (cu) and ($7), or [*lu) and 
Coulomb fields but NOT Coulomb force fields because they can NOT generate any direct 
local change of momentum, since as we mentioned earlier, these fields are static and 
they conserve their energy, momentum and angular momentum. The entire local force 
is given in the theory by the divergence of the Maxwell stress tensor which is a nonliner object, 
namely, a bilinear combination of the field components and their derivatives and all its terms 
are mutually compenseted in the static case. 

As we saw, the Coulomb force gets an admissible interpretation as an integral characteristic 
of the system describing some realizable integral tendences to minimization of the integral inter- 
action energy U of the two fields at the expense of the kinetic energy of the two particles. Surely, 
u! and carry some local physical information but in a quite indirect manner: except spherical 
symmetry (which, of course, is not specific only for electric fields) any of these two local objects 
can not clarifiy the physical nature of the local changes in the space when charged particles 
are around. In other words, from local point of view, we could not identify uj as electric field. 
Any topologically nontrivial region of the kind "M'^ minus a point" generates such fields, so, the 
electric nature of the field can be proved only by means of additional procedures concerning the 
integral structure of the system. 

The topological interpretation of to and 17 also suggests that the description is rather integral 
than local: although uo and ft are local objects, in fact they are just specially chosen representa- 
tives of integral characteristics of the physical system considered: they specify the topology of 
the space where the two fields are defined. For another example, the Newton gravitation force 
law looks the same except the different interpretation of the corresponding topological numbers 
as masses. Following the same argument, the Newton gravitation force law is of integral nature 
and shows similar tendences except that the masses are always positive numbers, so, the cor- 
responding interaction energy should be always negative, which does not allow repulsion. But 
this integral difference says too little about the local nature of the two physically different field 
structures. 

Remark. As for the relativistic formulation, uj and 0, can also be correctly defined and used, 
just the topology of the space where the two fields are defined very slightly changes: S S x M. 

The above consideration makes us think that, from theoretical point of view, the Maxwell 
stress tensor field is the right object in terms of which local force fields must be defined, namely, 
through computing its divergence. If the field is static and free then this divergence is zero 
and represents physically admissible quasistatic local changes, i.e. corresponding equations of 
motion, and any additional conditions must be consistent with this zero divergence. In our static 
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case we have da; = 0, d * a; = 0, so this divergence is zero: 



da;'' = rot E X E + E divE = 0, 



where ^lo^^ {dLo)ijkdx'' = EdivE, {*ioy{d * Lo)if;dx^ = (rotE) x E, {x^,x'^,x^) = {x,y,z), and 
vectors and covcctors are identified through the cuclidcan metric. Hence, we can not gain energy 
from the field E in local way. Therefore, the Coulomb force has not a local nature. 

The situation seriously changes when we are going to consider independent and self-consistent 
time-dependent and time-stable spatially finite physical objects of electromagnetic nature, namely, 
we have no such topologically motivated suggestions to choose adequate mathematical objects 
been able to represent appropriately the corresponding physical stresses. Hence, the mathemat- 
ical model must be created on the basis of assumptions of quite different nature, for example: 
requirements for definite and appropriately defined integrability properties representing the ob- 
ject's time stability; experimentally proved and traditionally assumed straight-line propagation 
of the energy-density; orthogonality of the electric and magnetic components of the field sug- 
gesting absence of local interaction energy between the electric and magnetic components or 
their new versions; some notion for internal energy redistribution during time-evolution, etc. In 
our view, in such cases, the eigen and other algebraic properties of the corresponding stress- 
energy-momentum tensor field should play a basic role. 

In view of this in the next two chapters we consider in a more detail from formal point of 
view the general Maxwell stress tensor as a starting physically meaningful theoretical object . 



2.2 A non-physical view on Maxwell stress tensor 

The mathematical identities have always attracted the attention of theorists, in particular, 
those identities which involve the derivatives of the objects of interest (differential identities). 
A well known such example is the Bianchi identity satisfied by any connection components: 
this identity is a second order system of (in general, nonlinear) partial differential equations. 
The gauge interpretation of classical Maxwell electrodynamics, as well as the Yang-Mills theory, 
substantially make use of this identity. Such identities are of particular importance when on the 
two sides of "=" stay correctly (i.e. in a coordinate free way) defined expressions. 

The physical reason to consider a couple of vector fields on M'^ (or a couple of spatial- 
directed vector fields on [M"^ x {time — coordinate)] as mathematical representation of a free time 
dependent physical field comes from the observation that, being free, it propagates translationaly 
along straight lines, so the identification properties of the field should be searched inside the 
2-dimensional distribution that is orthogonal to these straight lines. 

We begin with the well known differential relation satisfied by the square of every vector 
field V on the euclidean space M^. Our attention is directed to the square of V just because 
of the experimentally suggested assumption that mesures the energy-density of each of the 
electric and magnetic components. 

Let be related to the standard coordinates (x* = x,y,z),i = 1,2,3; we denote by " x " 
the vector product, and make use of the V-operator: 

lv(y2) = V xvotV + {V.V)V = VxTotV + VvV. 

Clearly, on the two sides of this relation stay well defined quantities, i.e. quantities defined 
in a coordinate free way. The first term on the right hand side of this identity accounts for 
the rotational component of the change of V, and the second term accounts mainly for the 
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translational component of the change of V. Making use of component notation we write down 
the last term on the right side as follows (summation over the repeated indices): 



Substituting into the first identity, and making some elementary transformations we obtain 



where 6"^^ = 1 for i = j, and (5*-^ = for i ^ j are the euclidean metric components. If now W is 
another vector field it must satisfy the same above identity: 



We note the invariance of M*^ with respect to the transformations (V, W) — > (Fcosa — 
Tysina, Vsina+Wcosa) where a{x, y, z) is an arbitrary real function. With respect to the slightly 
more general transformation (y, VF) {V a — W h,V h + W a) where (a, 6) are real nonzero 
functions, we obtain M{V,W) (a^ + b^)M{V,W). Hence, the transformations {V,W) 
{V a — W b,V b + W a) do not change the eigen directions structure of M*^. 

The expression inside the round brackets on the left of (2.3), denoted by M*-^, looks formally 
the same as the introduced by Maxwell tensor from physical considerations concerned with the 
electromagnetic stress energy properties of continuous media in presence of external electromag- 
netic field. This allows to call formally any such tensor Maxwell stress tensor generated 
by the two vector fields (V^,Ty). The term "stress" in this general mathematical setting is not 
topologically motivated as in the Coulomb case, but could be justified in the following way. 
Every vector field on generates corresponding flow by means of the trajectories started from 
some domain Uo C M^: at the moment t > the domain Uo is diffeomorphically transformed to 
a new domain Ut C M^. Having two vector fields on we obtain two compatible flows, so, the 
points of any domain Ug C are forced to accordingly move to new positions. 

We emphasize the following moments: first, the identity (2.3) is purely mathematical; sec- 
ond, on the two sides of (2.3) stay well defined coordinate free quantities; third, there is no 
V ^ W interaction stress: the full stress is a sum of the T^-stress and T^-stress. 

Physically, we say that the corresponding physical medium that occupies the spatial region 
Uo and is parametrized by the points of the mathematical subregion Uo C M^, is subject to 
compatible and admissible physical " stresses" generated by physical interactions mathematically 
described by the couple of vector fields (y, W), and these physical stresses are quantitatively 
described by the corresponding physical interpretation of the tensor M*-^(V, W). 

We note that the stress tensor M^^ in (2.3) is subject to the divergence operator, and if we 
interpret the components of M^^ as physical stresses, then the left hand side of (2.3) acquires 
in general the physical interpretation of force density. Of course, in the static situation as it 
is given by relation (2.3), no energy-momentum propagation is possible, so at every point the 
forces mutually compensate: VjM*-' = 0. If propagation is allowed then the force field is NOT 
zero: VjM*-' ^ 0, and we may identify the right hand side of (2.3) as a real time-change of 





Summing up these two identities we obtain the new identity 




[(rot V)xV + Vdw V + (rot W) x W + Wdiv Wy . 



(2.3) 
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appropriately defined momentum density P. So, assuming some expression for this momentum 
density P we are ready to write down corresponding field equation of motion of Newton type 
through equalizing the spatially directed force densities VjM*-^ with the momentum density 
changes along the time coordinate, i.e. equalizing VjM*-' with the rf-derivative of P, where 
c = const is the translational propagation velocity of the momentum density flow of the physical 
system (y, VF). In order to find how to choose P we have to turn to the intrinsic physical 
properties of the field, so, it seems natural to turn to the eigen properties of M*-^ , since, clearly, 
namely M^^ is assumed to carry the physical properties of the field. 



2.3 Notes on the eigen properties of Maxwell stress tensor 

We consider M*-'(E,B) at some point p G and assume that in general the vector fields E 
and B are linearly inependent, so E x B 7^ 0. Let the coordinate system be chosen such that 
the coordinate plane {x,y) to coincide with the plane defined by E(p),B(p). In this coordinate 
system E = (Ei,E2,0) and B = (i?i,i?2,0), so, identifying the contravariant and covariant 
indices through the Euclidean metric 5^^ (so that M*-^ = Mj = Mjj), we obtain the following 
nonzero components of the stress tensor: 

Ml = {E^f + {B^f - ^(E^ + B^); = Mf = E^E2 + Bi B^; 

M| = {E^f + (5^)2 - ^(E2 + B2); M| = -^(E2 + B2). 

Since = — M|, the trace of M is Tr{M) = — ^(E^ + B2). The eigen value equation acquires 
the simple form [(M^)^ — (A)^] + (Mg)^] (M3 — A) = 0. The corresponding eigen values are 



Ai = --(E2 + B2); A2,3 = ±V(^1 )' + (^2)' = ^l^^{h? + {h?, (2.4) 

where /i = B^ — E^, I2 = 2E.B. The corresponding to Ai eigen vector Zi must satisfy the 
equation E(E.Zi) + B(B.Zi) = 0, and since (E, B) are linearly independent, the two coefficients 
(E.Zi) and (B.Zi) must be equal to zero, therefore, Zi 7^ must be orthogonal to E and B, 
i.e. Zi must be colinear to E x B: 

The other two eigen vectors ^2,3 satisfy correspondingly the equations 

E(E.Z2,3) + B(B.Z2,3) = [ ± \^/{hY + {h? + ^(E2 + B2)] Z2,3. (*) 

Taking into account the easily verified relation 



{hf + {hf = [ ) - |E X B|^ so ^ |E X B| > , (2.5) 



E2 + b2\2 _ E^ + B 



1 



we conclude that the coefficient before ^2,3 on the right is always different from zero, therefore, 
the eigen vectors .^2,3 (p) he in the plane defined by (E(p), B(p)), p G M^. In particular, the above 
mentioned transformation properties of the Maxwell stress tensor M{V, W) {a? + b'^)M(y, W) 
show that the corresponding eigen directions do not change under the transformation (V, 1^) 
{Va-Wb,Vb + Wa). 

The above consideration suggests that the intrinsically determined potential dynamical abil- 
ities of propagation of the field are: translational along (E x B), and rotational inside the plane 
defined by (E, B). 
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It is natural to ask now under what conditions the very E and B may be eigen vectors of 
M(E, B)? Assuming A2 = (A)^ + (^2)^ and Z2 = E in the above relation (*) and having in 
view that E x B / we obtain that E(E2) + B(E.B) must be proportional to E, so, E.B = 0, 
i.e. I2 = 0. Moreover, substituting now /2 = in that same relation we obtain 

e2 = ^(B2-E2) + 1(e2 + b2) = b2, i.e., 7i = 0. 

The case "-" sign before the square root, i.e. A3 = —^^/(Jip~+Jl2p, leads to analogical 
conclusions just the role of E and B is exchanged. 

Corollary. E and B may be eigen vectors of M(E, B) only if Ii = I2 = 0. 

The above notices suggest to consider in a more detail the case A2 = — A3 = for the 
vacuum case. We shall show, making use of the Lorentz transformation in 3-dimensional form 
that, if these two relations do not hold then under E x B 7^ the translational velocity of 
propagation is less then the speed of light in vacuum c. Recall first the transformation laws of 
the electric and magnetic vectors under Lorentz transformation defined by the 3-velocity vector 

V and corresponding parameter /3 = v/c,v = \v\. If 7 denotes the factor 1/^/1-/3^ then we 
have 

E' = 7E + ^^v(E.v) + -V X B, 

c 

B' = 7B + ^-^v(B.v) - -v X E. 

c 

Assume first that I2 = 2E.B = 0, i.e. E and B are orthogonal, so, in general, in some 
coordinate system we shall have E x B 7^ . 

If /i > 0, i.e. |E| < |B|, we shall show that the conditions E' = 0,v.B = 0, 00 > 7 > 
are compatible. In fact, these assumptions lead to 7V.E + (1 — 7)(E.v) = 0, i.e. E.v = 0. 
Thus, c|E| = t;|B||sin(v,B)|, and since v.B = then |sin(v, B)| = 1. It follows that the speed 

IE 

V = cjgj < c is allowed. 

If Ii < 0, i.e. |E| > |B|, then the choice B' = and v.E = analogically lead to the 

IB 

conclusion that the speed v = cj^ < c is allowed. 

Assume now that I2 = 2E.B 7^ 0. We are looking for a reference frame K' such that 

E' X B' = 0, while in the reference frame K we have E x B 7^ 0. We choose the relative velocity 

V such that v.E = v.B = 0. Under these conditions the equation E' x B' = reduces to 

ExB + -(e2 + b2) =0, so, - = IE x BI/fE^ + b2). 

c c 

Now, from the above mentioned inequality E^ + B^ — 2|E x B| > it follows that ^ < 1. 

Physically, these considerations show that under nonzero Ii and I2 the translational velocity 
of propagation of the field, and of the field energy density of course, will be less than c. Hence, 
the only realistic choice for the vacuum case (where this velocity is assumed by definition to be 
equal to c), is Ii = I2 = 0, which is equivalent to E^ + B^ = 2|E x B|. Hence, assuming \Tr{M)\ 
to be the energy density of the field, the names "electromagnetic energy flux" for the quantity 
cE x B, and "momentum" for the quantity ^E x B, seem well justified without turning to any 
field equations. 

These considerations show also that if Ii = 0, i.e. |Ep = jBp during propagation, then 
the electric and magnetic components of the field carry always the same energy density, so, a 
local mutual energy exchange between E and B is not forbidden in general, but, if it takes 
place, it must be simultanious and in equal quantities. Hence, under zero invariants Ii = 
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and I2 = 2E.B = 0, internal energy redistribution among possible subsystems of the field is 
allowed but such an exchange should occur without available interaction energy because the full 
energy density is always equal to the sum of the energy carried by the electric and magnetic 
components of the field. However, the required time stability and propagation with velocity 
"c" of the field suggest/imply also available internal momentum exchange since under these 
conditions the energy density is always equal to the momentum magnitude |E x B|, and E and 
B can not carry momentum separately. Moreover, besides (E, B), another subsystem of the 
field has to be constructed out of (E, B) such, that both these two subsystems to carry always 
the same quantity of energy-momentum, so the exchange also must be in equal quantities. 

After these preliminary considerations we procedc to write down dynamical equations for the 
field through specializing how the internal local momentum exchange is realized keeping always 
in mind that the free field energy density propagates translationaly with the speed of light, so 
the relations 

h=l2 = 0, i.e. E^ + B^ = 2|E x B| 

must always hold. 

2.4 Nonlinear equations for the electromagnetic field 

We are going to consider time dependent fields, and begin with noting once again that the 
assumption that the energy density of the field coincides with |Tr(M)| = ^[E^ + B^] presupposes 
that there is NO interaction energy between the electric and magnetic components of the field: 
the full stress tensor (and the energy density, in particular) is a sum of the stress tensors 
determined separately by E and B. Of course, this does NOT mean that there is no energy 
exchange between the electric and magnetic components, but if such an exchange takes place, 
it must occur simultaniously and in equal quantities.. 

Now, following the above stated idea that the field momentum density could be responsi- 
ble for such an internal energy-momentum exchange, we have to find at least two appropriate 
subsystems of the field which subsystems are NOT the electric E and magnetic B ones, but 
are constructed out of them. Note that the assumption that the field momentum is given by 
^E X B, i.e. it is a bilinear function of the electric and magnetic components, and that the 
local energy is always equal to |E x B|, suggests that the electromagnetic momentum of the 
field is of interaction nature. The point now is to get some clarification how such a local mo- 
mentum exchange (and the corresponding energy exchange) takes place and to find appropriate 
mathematical representatives of the corresponding partners realizing such special kind of energy 
exchange, since neither E nor B are able to carry momentum separately (although each of them 
may carry energy independently of the other). 

In view of the above we shall assume that the field keeps its identity through adopting some 
special and appropriate dynamical behavior according to its intrinsic capabilities. Hence, the 
corresponding dynamical/field equations must be consistent with the intrinsic stress-energy- 
momcntum nature of the field. So, our basic assumption is that the Maxwell stress tensor 
M(E, B) should play the basic role, and its zero- divergence in the static case should suggest 
how to determine the allowed dynamics. 

Recalling that any member of the family 

(£^,i3) = (Ecosa — Bsina; E sin a -|- B cos a) , a = a{x,y, z,t), 

generate the same Maxwell stress tensor, the most natural assumption should read like this: 
the field (E, B) is looking for an energy-momentum exchanging partner inside the 
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a{x, y, z, t)-familly of (E, B)-couples, and any such couple identifies itself through ap- 
propriate (local) interaction, defining in this way corresponding dynamical behavior. 

Replacing (V, W) in (2.3) with (E, B) we obtain 

/ E^ + B^' 

ViM'^ = Vi E'E^' + B^B^' - S'^ - 



2 

= [(rot E) X E + Ediv E + (rot B) x B + Bdiv B] ^ (2.6) 

As we mentioned, in the static case, i.e. when the vector fields (E, B) do not depend on the 
time coordinate ^ = ct, NO propagation of field momentum density P should take place, so, 
at every point, where (E, B) / 0, the stress generated forces must mutually compensate, i.e. 
the divergence VjM*-' should be equal to zero: VjM*'' = 0. In this static case Maxwell vacuum 
equations 

5B (?E 
rotE+ — =0, rotB- — =0, divE = 0, divB = (*) 

give: rotE = rotB = 0; divE = divB — 0, so, all static solutions to Maxwell equations determine 
a sufficient, but NOT necessary, condition that brings to zero the right hand side of (2.6) through 
forcing each of the four vectors there to get zero values. 

In the non-static case, i.e. when ^ / 0; ^ 7^ 0, time change and propagation of field 
momentum density should take place, so, a full mutual compensation of the generated by the 
Maxwell stresses at every spatial point local forces may NOT be possible, which means ViM'-' 7^ 
in general. These local forces generate time-dependent momentum propagation at the spatial 
points. Therefore, if we want to describe this physical process of field energy-momentum density 
time change and spatial propagation we have to introduce explicitly the dependence P(E,B). 
If we follow the classical (nonrelativistic) way of consideration and denote by 5 the vector field 
with components = ViM^^ , we can write down the force flow across some finite 2-surface S 
in the usual (and widely spread in almost all textbooks) way as Jg ^.ds (from modern point 
of view we should write i:g{dx A dy A dz) instead of 5-ds under the integral, where denotes 
the inner product between the vector field 5 and the volume form dx A dy A dz, i.e. to make 
use of the Poincare isomorphism between vector fields and 2- forms on R^). This flow generates 
changes of the momentum density flow across S which should be equal to ^ Jg P(E, B).ds. We 
obtain 

^ / P(E,B).ds= / 5.ds. 
dt Js Js 

The explicit expression for P(E, B), paying due respect to J.Poynting [13], and to J.J.Thomson, 
H. Poincare, M. Abraham [15], and in view of the huge, a century and a half available experence, 
has to be introduced by the following 

Assumption: The entire field momentum density is given byP :=IExB . 

According to the Assumption and the above interpretation of the relation ViAP^ 7^ 0, and 
in view of the arbitrariness of the 2-surface S we come to the vector differential equation 

^(ExB)=^, C = ct, (**) 

which according to relation (2.6) is equivalent to 

9B \ / 9E \ 

Tot'E + — \ X E-KEdivE-F frotB - — j x B-^BdivB = 0. (2.7) 
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This last equation (2.7) we write down in the following equivalent way: 



dB , 

rotE + — 1 X E + BdivB = - 
9^ 



5E 

rotB- — 1 X B + EdivE 
9^ 



(2.8) 



The above relation (**) and the corresponding differential relation (2.7)/(2.8) we consider as 

mathematical adequate in momentum-change terms of the electric-magnetic and magnetic- 
electric induction phenomena in the charge free case. We recall that these induction phenomena 
are described in what we call "Faraday-Maxwell theory" by the following well known integral 
and differential equations 



^ / B.ds = - / rotJii.ds ^ 

Js Js 9^ 

d_ 

JS JS 

We stress once again that these last Faraday-Maxwell relations have NO direct energy-momentum 
change-propagation (i.e. force flow) nature, so they could not be experimentally verified in a 
direct way. Our feeling is that, in fact, they are stronger than needed. So, on the corresponding 
solutions of these equations we'll be able to write down formally adequate energy-momentum 
change expressions, but the correspondence of these expressions with the experiment will cru- 
cially depend on the nature of these solutions. As we already mentioned, the nature of the free 
solutions (with no boundary conditions) to Maxwell vacuum equations with spatially finite and 
smooth enough initial conditions requires strong time-instability (the Poisson theorem for the 
D'Alembert wave equation). And time-stability of time-dependent vacuum sohitions Tisually 
requires spatial infinity (plane waves), which is physically senseless. Making calculations with 
spatially finite parts of these spatially infinite solutions may be practically acceptable, but from 
theoretical viewpoint assuming these equations for basic ones seems not acceptable since the 
relation "time stable physical object - exact free solution" is strongly violated. 

Before to go ftirther we write down the right hand side bracket expression of (2.8) in the 
following two equivalent ways: 



rotE.ds 



dB 



= —rotE, (the Faraday induction law). 



rotB.ds 



^=rotB, 



(the Maxwell displacement current law). 



rotB 



X B + (-E)div (-E) 



9E \ 

rot (-B) + — j X (-B) + EdivE 



(2.9) 



These last two expressions (2.9) can be considered as obtained from the left hand side of (2.8) 
under the substitutions (E, B) (B, — E) and (E, B) (— B,E) respectively. Hence, the field 
(E, B) chooses as a partner-field one of the fields (— B, E), or (B, — E). 
We may resume this in the following way: 

An adequate mathematical representation of a time dependent free electromag- 



(E,B);(-B,E) 



or 



(E,B);(B,-E) 



netic field requires a collection of two fields : 

We could also say that a real free field consists of tw^o interacting subsystems described 
by two partner-fields inside the a{x, y, z, i)-family 

{S,B) = (E cos a — B sin a; E sin a -|- B cos a) 

giving the same Maxwell stress-energy tensor. Each partner-field has electric and 

magnetic components, and each partner-field is determined by the other through 
(±^) - rotation-like transformation. Both partner-fields carry the same energy- 
momentum and minimize the relation -|- /| ^ 0. This view and relation (2.7/2.8) sug- 
gest, in turn, that the intrinsic dynamics of free real time-dependent electromagnetic 
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fields could be considered as establishing and maintaining local energy-momentum 
exchange partnership between two fields called above partner-fields, and, since (E, B) 
and (— B, E) carry always the same stress-cncrgy-momentum, the allowed inter-exchange 
is necessarily simultaneous and in equal quantities, so, each partner-field conserves 
its energy-momentum. 

In order to find how much is the locally exchanged energy-momentum we are going to inter- 
pret the equation (2.8) in accordance with the view on equations of motion as stated in Sec. 2. 2. 

Our object of interest representing the integrity of a real time dependent electromagnetic 



field, is the couple 



fE,B) 



B,E1 



fE,B);(B,-El 



is considered analogi- 



(the other case 

cally). In view of the above considerations our equations should directly describe admissible 
energy-momentum exchange between these recognized two subsystems, i.e. from formal point 
of view, between the two partner-fields. Hence, we have to define the corresponding change- 
objects D(E, B) and D(— B,E) for each partner-field, and their self-" projections" and their 
mutual "projections". 

The change object -D(E, B) for the first partner-field (E, B) we naturally define as 

(SB \ 
rotE + -^; divBj . 

The corresponding "projection" of I?(E, B) on (E, B) 

<P[I)(E,B);(E,B)] =^ 

is suggested by the left hand side of (8) and we define it by : 

" / 5B \ 1 / 5B \ 

<P (rotE + — ; divB j ; (E, B) := (rot E + —j x E Bdiv B. 

For the second partner-field (— B,E), following the same procedure we obtain: 



(SB \ 
rotE+-^; divBj ; (E, B) 



^[Z?(-B,E);(-B,E)]=q} 

dE 



dE 

rot(-B) + — ; divE) ;(-B,E) 



rot (-B) + 



X (-B) -FEdivE = rotB 



dE 



X B -I- Ediv E. 



Hence, relation (2.7) looks like 

[DiE, B); (E, B)] + ^ [D{-B, E); (-B, E)] = 0. 

The accepted two-component view on a real time dependent electromagnetic field allows 
in principle admissible energy-momentum exchange with the outside world through any of the 
two partner-fields. Hence, the above calculations suggest to interpret the two sides of (2.8) 
as momentum quantities that each partner-field (E, B), or (— B,E), is potentially able to give 
to some other physical object withot destroying itself, and these quantities are expressed in 
terms of E, B and their derivatives only. In the case of free field, since no energy-momentum 
is lost by the field, there are two possibilities: first, there is NO energy-momentum exchange 
between the two partner-fields, second, each of the partner-fields changes its energy-momentum 
at the expense of the other through simultanious and in equal quantities exchanges. Such 
kind of mutual exchange is in correspondence with the mathematical representatives of the 
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two subsystems: the partner-fields (E, B) and (— B, E) being members of the above mentioned 
a{x, y, z, f)-family, obviously carry the same energy |(E^ + B^) and momentum ^(E x B). If we 
denote by An and by A22 the allowed energy-momentum changes of the two component-fields, 
by A12 the energy-momentum that the first partner-field receives from the second partner-field, 
and by A21 the energy-momentum that the second partner-field receives from the first partner- 
field, then according to the energy-momentum local conservation law we may write the following 
equations: 

Aii = Ai2 + A2i; A22 = - (A21 + A12) , 

which is in accordance with the equation (2.8): An + A22 = 0. 

We determine now how the mutual momentum exchange between the two partner-fields 
P(E,B) P(-B,E)) 01') -^(E.B) ^ P(B,-E) is performed, i.e. the explicit expressions for A12 and 
A21, keeping in mind that |P(e,b)I = |P(-B,E)I = |P(B-E)I = I^P[{E,B);(-B,E)]I- The formal 
expressions are easy to obtain. In fact, in the case P(e,b) P(-b,e)) i-e- the quantity A21, we 
have to "project" the change object for the second partner-field given by 



Z?(-B,E) := (^rot(-B) + ^; divE^ 



on the first partner-field (E, B). We obtain: 

/ 9E\ / 0E\ 

A21 = frot(-B) + — j X E-F BdivE = - frotB-— j xE-F BdivE . (2.10) 

In the reverse case P(_b,e) ~^ P(E,B)i i-^- the quantity A12, we have to project the change-object 
for the first partner-field (E, B) given by 

Z)(E,B) := ^rotE+ divB^ 

on the second partner-field (— B,E). We obtain 

(9B \ / 9B \ 

rot E + — J x (-B) + Ediv B = - (rot E + — \ x B Ediv B. (2.11) 

So, the internal local momentum balance is governed by the equations 



' (9B , 

rotE-l- — 1 X E-FBdivB = 



- rot E 



5B\ „ „ / „ 9E 

— jxB + EdlvB-(^rotB- — 



xE-|-BdivE, 
(2.12) 



(5E \ / 9E \ / c?B \ 

rotB - — J X B-l-EdivE = (rotB - — j x E- BdivE -|- f rotE-F — J x B - EdivB. 

(2.13) 

These two vector equations (2.12)-(2.13) we consider as natural Newton type field equations. 
According to them the intrinsic dynamics of a free electromagnetic field is described by two 
couples of vector fields, [(E, B); (— B, E)], or [(E, B); (B, — E)], and this intrinsic dynamics could 
be interpreted as a direct energy-momentum exchange between two appropriately individualized 
subsystems mathematically described by these two partner-fields. 
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A further natural specilization of the above two vector equations (2. 12)- (2. 13) could be made 
if we recall that this internal energy-momentum exchange realizes a special kind of dynamical 
equilibrium between the two partner-fields, namely, the two partner-fields necessarily carry 
always the same energy and momentum : M*-'(E, B) = M^^{£,B), so each partner-field must 
conserve its momentum : An = A22 = 0. In such a dynamical situation each partner-field loses 
as much as it gains during any time period, so, equations (2.12)-(2.13) reduce to 

Aii= rotE-F— xE-FBdivB = 0, (2.14) 



A22 = (^rot B - X B + Ediv E = 0, (2.15) 

3B \ / 9E \ 

A12 + A21 = ( rotE-l- — I X B - EdivB -I- ( rotB - — ) x E - BdivE = 0. (2.16) 



Equation (2.16) fixes, namely, that the exchange of energy-momentum density between the two 
partner-fields is simultanious and in equal quantities, i.e. a permanent dynamical equi- 
librium between the two partner-fields holds: P(e,b) P(-b,e)) or, P(e,b) ^ P(b,-E)- 

Note that, if equations (2.14) and (2.15) may be considered as field-equivalents to the zero 
force field (eqn. (2.14)) and its dual (eqn. (2.15)), this double-field viewpoint and the corre- 
sponding mutual energy-momentum exchange described by equation (2.16) are essentially new 
moments. Equations (2.14)-(2.16) also suggest that the corresponding fields are able to exchange 
energy-momentum with other physical systems in three ways. If such an exchange has been done, 
then the exchanged energy-momentum quantities can be given in terms of the characteristics 
of the other physical system (or in terms of the characteristics of the both systems) and to be 
correspondingly equalized to the left hand sides of equations (2.14)-(2.16) in accordance with 
the local energy-momentum conservation law. 

2.4.1 Some Properties of the nonlinectr solutions 

Clearly, all solutions to Maxwell pure field equations (*) are solutions to our nonlinear equations 
(2.14)-(2.16), we shall call these solutions linear, and will not be interested of them just because 
the notion for stress-energy-momentum partnership between (E,B) and (— B,E) is missing. 
Therefore, we shall concentrate on those solutions of (2.14)-(2.16) which satisfy the conditions 

(9B (?E 
rotE-F — 7^0, rotB- — 7^0, divET^O, divBT^O. 

These solutions we call further nonlinear (among them there are no constant ones as it is in the 
class of linear ones). We note some of the properties they have. 

1. E.B = 0; 

2. (^rotE + .B = 0; (^rot B - ^) .E = 0. 

Prom these two relations the classical Poynting energy-momentum balance equation follows. 

The above two properties are obvious from equations (2.14) and (2.15). 

3. If (E, B) defines a solution then (E', B') = (aE - 6B; 6E -|- aB), where a, 6 e M, defines 
also a solution. This property is immediately verified through substitution. 

4. E^ = B^. 
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To prove this, we first multiply equation (2.14) on the left by E and equation (2.16) by B (scalar 
products). Then we make use of the above properties 1 and 2, of the vector algebra relation 
X.{Y X Z) = Z.(X X Y), and of the assumed nonlinear values of the divergences of E and B. 

Properties 1. and 4. say that all nonlinear solutions to (2.14)-(2.16) are null fields, i.e. the 
two well known relativistic invariants Ji = B^ — E^ and I2 = 2E.B of the field are zero, and 
this property leads to optimisation of the inequality if + I2 ^0 (recall the cigen properties of 
Maxwell stress tensor. Sec. 2. 3), which in turn guarantees a(a:;, y, z, i)-invariance of Ii and l2- 

5. B. (rot B - - E. (rot E + ^) = B.rotB - E.rotE = 0. 

To prove this property we first multiply (vector product) (2.14) from the right by E, recall 
property 1., then multiply (scalar product) from the left by E, recall again E.B = 0, then 
multiply from the right (scalar product) by B and recall property 4. 

Property 5. suggests the following consideration. If V is an arbitrary vector field on then 
the quantity V.rotV is known as local helicity and its integral over the whole (compact) region 
occupied by V is known as integral helicity, or just as helicity of V. Hence, property 5. says 
that the electric and magnetic components of a nonlinear solution generate the same helicities. 
If we consider (through the euclidean metric) E as 1-form on and denote by d the exterior 
derivative on , then E A dE = E.rotE dx Ady A dz, so, the zero helicity says that the 1-form 
E defines a completely integrable Pfaff system: E A dE = 0. The nonzero helicity says that the 
1-form E defines non-integrable Id Pfaff system, so the nonzero helicity defines corresponding 
curvature. Therefore the equality between the E-helicity and the B-helicity suggests to consider 
the corresponding integral helicities Jjg3 E A dE = Jj^g B A dB (when they take finite nonzero 
values) as a measure of the spin properties of the solution. 



6. Example of nonlinear solution: 
E 



z z 
(p{x^ ct + ez)cos{~K,- — h const)^ ^(x, ct + ez)sm{—K,— + const)^ 



B = 



z z 
£(f){x, y, ct + £z) sm{—K- — h const), —e(f){x, y,ct + £z)cos{—k- — h const), 



where cj){x, y, ct + ez) is an arbitrary positive function, < oo is an arbitrary positive constant 
with physical dimension of length, and £ and k take values ±1 independently. Modifying the 
helicity 3-forms to k-^E A dE = k-^B A dB, then the corresponding 3d integral gives kTE, 
where k = ±1, T = Al^/c and E = J (fp'dxdydz is the integral energy of the solution. 



2.5 Discussion 

The main idea of this part of the paper is that carrying out the Newton way for writing down 
dynamical equations for particles in mechanics to writing down dynamical equations for con- 
tinuous field systems should naturally result to nonlinear partial differential equations even in 
non-relativistic theories. Moreover, clarifying the sense of the information included in these 
dynamical equations according to the Newton approach, we come to the conclusion formulated 
in the Introduction, namely, we have to mathematically describe those changes of the object 
considered which are qualified as admissible and consistent with the system's identification and 
with the local energy- momentum balance relations. In the case of "free" systems these re- 
lations represent the local energy-momentum exchange/conservation properties of the system. 
The energy-momentum characteristics are chosen because of their two important properties: 
they are physically universal and conservative. This means that every physical object carries 
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nonzero energy-momentum and, vice versa, every quantity of energy-momentum is carried by 
some physical object. Hence, if a physical object loses/gains some quantity of energy-momentum 
then some other physical object necessarily gains/loses the same quantity of energy- momentum. 
If this viewpoint is assumed, then the problem of finding appropriate dynamical equations for an 
object reduces mainly to: first, getting knowledge of the potential abilities of the object consid- 
ered to lose and gain energy-momentum; second, to create adequate mathematical quantities 
describing locally these abilities. 

The electromagnetic field, considered as a continuous physical object of special kind, gives a 
good example in this direction since, thanks to Maxwell's fundamental and summarizing works, 
all the information needed is available. The notices of Poynting [13], and Thomson, Poincare 
and Abraham [14], showing the importance of the (deduced from Maxwell equations) vector 
-El X B from local energy-momentum propagation point of view, has completed the resource 
of adequate and appropriate mathematical objects since it appears as natural complement of 
Maxwell stress tensor, and allows to write down dynamical field equations having direct local 
energy-momentum balance sense. However, looking back in time, we see that this viewpoint for 
writing down field equations has been neglected, theorists have paid more respect and attention 
to the "linear part" of Maxwell theory, enjoying, for example, the exact but not realistic, and 
even physically senseless in many respects, plane wave solutions in the pure field case. 

Therefore, not so long after the appearance of Maxwell equations the photoeffect exper- 
iments showed the nonadequateness of the linear part of Maxwell theory as a mathematical 
model of electromagnetic fields producing realistic model-solutions of free time-dependent fields. 
Although the almost a century long time development of standard quantum and relativistic 
quantum theories that followed, a reasonable model-solutions describing individual photons, 
considered as basic, spatially finite and time-stable objects, these theories have not presented so 
far. Nobody doubts nowadays that photons really exist, and this very fact suggests to try first 
classical field approach in finding equations admitting 3d-finite and time stable solutions with 
appropriate properties. 

The historical perspective suggests to follow the 4-potential approach, but modern knowledge 
and experience, and even the Maxwell stress tensor achievements, suggest some different views. 
In fact, we have all reasons to consider the microobjects as real as all other physical objects, so, no 
point-like charges and infinite field model-solutions should be considered as adequate. Since the 
4-potential approach in Maxwell theory docs not allow spatially finite and time stable pure field 
solutions with photon-like structure and behavior its interpretation as a basic concept does not 
seem to be appreciable. Also, the 4-potential approach excludes many solutions of the charge free 
Maxwell equations. For example, in relativistic terms the well known field F = ^dr Ad,^, dF = 
0, has global 4-potential, and its Minkowski-dual *F = qsin 9d9 Adip, d*F = 0, has NO global 
4-potential. Now, the 2-parameter family of 2-forms = {aF — b* F;bF + a* F),a,b eR, 

gives an infinite number of solutions to Maxwell equations = 0, d * 3^ = admitting NO 
global 4-potential. This suggests the view that the 4-potential can be used as a working tool 
(wherever it causes no controversies) but not as a basic concept. 

In conclusion, paying due respect to the Newton view on dynamical equations and to the 
local energy-momentum conservation law we based our approach on the Maxwell stress tensor 
and on the Poynting vector as natural quantities carrying the physically meaningful energy- 
momentum characteristics of the electromagnetic field. The natural description in these terms 
leads to the assumption that any real time-dependent electromagnetic field consists of two in- 
teracting subsystems mathematically represented by the two partner- fields: [(E,B)], [(— B,E)]/ 
[(E, B)], [(B, — E)], or any couple [{£,B),{—B,£)] inside the considered a(x, y, z, t)-family of 
fields. These two partner-fields carry always the same stress-energy-momentum, and a dy- 
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namical equilibrium between these two subsystems is realized through a simultanious mutual 
energy-momentum exchange in equal quantities. The equations obtained represent formally this 
dynamical equilibrium, i.e. they show that partner-fields identify/recognize each other 
through appropriate local energy-momentum exchange partnership minimizing the 
quantity + I2 = ^(E^ + B^) — |E x B| ^0. Accordinly, all nonlinear solutions have zero 

invariants Ii = I2 = 0, and can not be constant. Among these zero-invariant nonlinear solutions 
there are time-stable and spatially finite ones with helical spatial structure, having photon-like 
properties and behavior. An analog of the Planck relation E = hv holds for these solutions, 
where the constant h appears as an integral helicity of such a solution. 
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Chapter 3 

Relativistic Considerations 



3.1 The Extended Electrodynamics approach 

The generalization of Classical Electrodynamics (CED) known as Extended Electrodynamics 
(EED) [11], starts with the conviction that CED surely carries inside the potential ability to be 
extended in such a way, that spatially finite and time-stable solutions of photon-like nature to 
be incorporated, and it exploits mainly two ideas: the idea for a direct local energy-momentum 
exchange sense of the new dynamical equations, and the well known dual symmetry (mentioned 
above) of the vacuum CED-equations and local conservation laws in the frame of relativistic 
formalism. Let's consider first the elementary physical approach. 
Maxwell vacuum equations 

dB dE 
rotE-F — =0, rotB- — = 0, divE = 0, divB = (*) 

clearly suggest that a free electromagnetic field has two vector components: electric E and 
magnetic B. On one hand, from physical viewpoint, these equations imply also interaction, 
i.e. energy-momentum exchange between the electric and magnetic components of the electro- 
magnetic field, and on the other hand, the energy density expression ^(E^ -|- B^) in the theory 
does not contain interaction energy term: the full energy density is the sum of the electric and 
magnetic energy densities. Now, since any of these two components can NOT carry momentum 
separately (the field momentum is given by ^E x B), then E and B can exchange only energy and 
NO momentum. But the field propagates translationally along null straight-line directions with 
the speed of light " c" , so it necessarily carries momentum being numerically equal (in energy 
units) to the energy-density, the so called electromagnetic energy flux. Hence, the field energy is 
of entirely dynamical nature and any internal energy exchange between subsystems necessarily 
implies corresponding momentum exchange. Therefore, from energy-momentum exchange point 
of view, we should be interested in finding such subsystems of the field, which are able to carry 
and exchange simultaniously both energy and momentum. 

In order to come to appropriate mathematical representatives in relativistic terms of such 
subsystems we recall that under null character of the local energy-momentum, i.e. when 
Tjj,yT^^ = 0, the translational propagation requires zero invariants: /i = B^ — E^ = and 
I2 = 2E.B = 0, which in terms of the relativistic 2-form formalism is equivalent respectively 
to h = \Fo,pF°'^ = and h = ^Fa(s{*F)'^^ = 0, where * is defined by the Minkowski 
(pscudo)metric: aA*(3 = —r]{a, P) y^]det(r]^dx^ A dx^ A dx^ A dx^. Now, recall the well known 
identity, being in force for any two 2-forms F and G in Minkowski space-time (M, ry) : 



iF^pG'^f'S-;, = F^.G^'^ - {*G)^,{*Fr. 
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Under G = F and ^F^/jF"^ = 0, it follows F^^-F'"" = (*F)^„(*F)'^''. In view of the canonical 
stress-energy-momentum tensor of the field (we omit the coefficient ^) 



rpu \ 

n - 2 



^F^pF-^S';^ - F^F'^- = -UF)^p{*FrHl - {*F),,{*FY-, 



4 M f*" 4' 

satisfying the Rainich condition T^T""" = {Taf^T"''^ S^^ = + I2W^{= in our case), this 
physically means that F(E, B) and *F(— B, E) can carry both energy and momentum, moreover, 
which is very important, they carry always the same stress-energy-momentum. Formally, this is 
partially hidden in the obvious invariance of with respect to F ^ *F. Now, since there is no 
other physical object participating in the energy-momentum exchange, we come to the conclusion 
that for adequate mathematical representatives of the two subsystems we are looking for, namely 
F and *F can be chosen. Moreover, since they necessarily carry always the same stress-energy- 
momentum, they may exchange locally energy-momentum only simultaniously and in equal 
quantities. Therefore, in view of V,.r^ = F"'^(dF)„^^ (*F)"'^(d * F)^/^,, = 0, a < P, the 
most natural dynamical equations should read 

F-^{dF)^p^ = 0, (*Fr^(d*F),;3^ = 0, (*Fr^(dF),;3^ + F"^(d*F),;3^ = 0, a < (3, 

where the first two equations require that F and *F conserve the energy-momentum they 
carry, and the third equation establishes the local dynamical equilibrium between F and *F: 
(>i:_p)"/^(dF)Q,^^ and F°'^(d * F)ad^ denote respectively the allowed from the consevation laws 
energy-momentum gains and losses of *F and F, which gains and losses are forbidden by the 
old equations dF = 0, d * F = 0. 

Extended Electrodynamics gives the following mathematical picture of this field dynamics. 
Recall that if (F, *F) is a CED vacuum solution, i.e. dF = 0, d * F = 0, then the combinations 
J^ = aF — b*F, T* = bF + a*F, where (a, 5) are two arbitrary real numbers, also give a CED 
vacuum solution and, since on Minkowski space the reduced to 2-forms Hodge star * satisfies the 
relation *^ = — i(iA2(Af)7 'wc obtain J^* = *T . The two corresponding energy tensors are related 
by T(T ^ T*) = {a^ + h"^) T{F, *F). Recall the real representation of complex numbers z = aI-\-bJ 
where I is the unit matrix in and J is the standard complex structure matrix in with 
columns (0, —1); (1,0). So, we obtain an action of the linear group of matrices a = al -\- bJ 
on the CED vacuum solutions. This is a commutative group G and its Lie algebra Q just adds 
the zero (2 x 2) matrix to G, and (/, J) define a natural basis of Q. So, having a CED vacuum 
solution, we have in fact a 2-parameter family of vacuum solutions. Hence, we can define a 
^-valued 2-form J7 on M by = F (g) / -|- *F J, and the equation dO = is equivalent to 
dF = 0, d * F = 0. 

Consider the new basis (/', J') oi Q given by 

l' = {aI + bJ), J' = {-bI + aJ). 

Accordingly, the "new" solution i.e. the old solution in the new basis of Q, will be 

n' = F(8)/'-F*F(8) J' = F®{aI-^bJ)-\-*F®{-bI-^aJ) = {aF-b*F)^I-\-{bF-^a*F)^J. 

In view of this we may consider this transformation as nonessential, i.e. we may consider (F, *F) 
and (JT, JT*) as two different representations in corresponding bases of Q of the same solution. 

Such an interpretation is approporiate and useful if the field shows some invariant properties 
with respect to this class of transformations. For example, if the Lorentz invariants 

h = l F^.Fi^'^ = (B2 - e2), 72 = 1 F^^i^Fr'' = 2E.B, 
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are zero: I\ = I2 = 0, then all the above transformations keep unchanged these zero-values of 
Ii and l2- In fact, under such a transformation [F,*F) — > {J-,J-*) the two Lorentz invariants 
transform to (/(j/g) in the following way: 

I[ = (a^ - + 2abl2, I2 = -2ahh + (a^ - b^) h, 

and the determinant of this transformation is (a^ + 6^)2 _^ q g^^ ^^^jj ggj^j stays a null field 
under these transformations. Moreover, NO non-null field can be transformed to a null field 
by means of these transformations, and, conversely, NO null field can be transformed to a non- 
null field in this way. Hence, the Lorentz invariance and the dual invariance of Ii and I2 hold 
simultanoiusly only in the null-field case. This observation distinguishes once again the null-field 
case. 

In order to come to the new equations we recall that every bilinear map ip : Q x Q ^ W, 
where W is some linear space with basis {ei},i = 1,2, ... , defines corresponding product in the 
^-valued differential forms by means of the relation 

(p{0.\ (g) Ci, (g) Cj) = Q\ A (g) ip{ei,ej). 

Now, let if = \/, where "V" is the symmetrized tensor product in Q. We consider the expression 

y{VL,*d9.). 

V(ri, *dJ^) = (F A *dF) (g)/V/-|-(*FA*d*F)(g)JVJ-F(FA*d*F(8)-F*FA *dF) ®I\J J. 
The vacuum EED equations are V(r2, *dr2) = 0, or equavalently, 

FA*dF = 0, (*F)A*d*F = 0, F A *d * F A *dF = 0. 
In terms of the codifferential 5 = *d* these equations look like 

6*F AF = 0, 5F A*F = 0, dF A F - 6 * F A *F = 0. 
In components we obtain correspondingly 

F"^(dF)„^^^(*FV(5*Fr = 0, (*Fr^(d*F)„^^^F^,(5F)^ = 0, a < (3; 

(*F)°'3(dF)„^^ + F"/5(d * F),^^ = (,5 * FYF,^ + {6Fn*F),^ = 0, a < p. 

It is easy to see that these equations are equivalent to the equations An = A22 = 0, A12 + A21 = 
as given in terms of (E, B) in Sec. 2. 4. Moreover, all nonlinear solutions to these EED vacuum 
equations, i.e. those satisfying dF 7^ 0, d * F 7^ 0, have zero invariants: /i = I2 = (for the 
case EED in presence of media see [11]). 

As for the energy-momentum tensor T^jy of the vacuum solutions, considered as a symmetric 
2-form on M, it is defined in terms of O as follows: 



T{X,Y) = ^*g[i{X)n,*i{Y)n] = -^X'^y'^[F^^F,^ + (*FV(*F),^' 



X^Y'^T^,, 



where {X, Y) are two arbitrary vector fields on M, g is the metric in Q defined by g{a, (3) = 
^tr{a.p*), and /?* is the transposed to (3. Note that g{I, J) = g{al + bJ, —hi — aJ) = 0, which 
elliminates the corresponding coefficient in T{X,Y), which reads F^(j{*FY'^ + {■^F)^„F'^'^ = 
5Fq,^(*F)"^(5^, so, in a (/-NONorthogonal basis of Q this coefficient will appear. Now, since 
^Fai3{*F)°'^ = 2(E.B), choosing ^-orthogonal basis in Q corresponds to mutual orthogonality 
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of (E, B) in this context, and choosing 5-nonorthogonal basis of G will formally result in some 
interaction between E and B. 

Finally, recall the generalization of Lie derivative Ck with respect to the fc-vector K, acting 
in the exterior algebra of differential forms according to the formula Ck = i{K)d — (— l)^di(K) 
[12]. Then, in view of the relations F^^F^*^ = (*F)^^F'*'^ = 0, the above equations acquire the 
form 

CpF = 0, £,-f{*F)=0, £j?(*F) + £(,~^)F = 0, 

where F and *F arc the ry-corresponding 2-vectors. In terms of ^} and Cl = F i^i ei + *F (8) 62 
these three equations can be represented as one relation as follows: 

= £pF ^ eiW ei + * F ^ €2 + {Cp * F + C^pF) ® ei V 62 = 0. 

The above consideration is based on the assumption that the ^-valued 2-form ri(F, *F) 
represents mathematically the wholeness+structural integrity of the electromagnetic field 
through the equations = 0, and direct physical motivation for such an assumption was not 
given. In order to motivate looking for photon-like solutions of these equations we now present 
our notion for photon-like object (s), and further we show how this notion leads to corresponding 
mathematics by means of which we could separate the desired subclass of solutions. 

3.2 The Notion of Photon-like Object(s) 
3.2.1 Introduction 

At the very dawn of the 20th century Planck [16] proposed and a little bit later Einstein [17] 
appropriately used the well known and widely used through the whole last century simple formula 
E = hiy, h = const > 0. This formula marked the beginning of a new era and became a 
real symbol of the physical science during the following years. According to the Einstein's 
interpretation it gives the full energy E of really existing light quanta of frequency u = const, 
and in this way a new understanding of the nature of the electromagnetic field was introduced: 
the field has structure, which contradicts the description given by Maxwell vacuum equations. 
After Dc Broglie's [18] suggestion for the particle-wave nature of the electron obeying the same 
energy-frequency relation, one could read Planck's formula in the following way: there are 
physical objects in Nature the very existence of which is strongly connected to some periodic 
(with time period T = l/i') process of intrinsic for the object nature and such that the Lorentz 
invariant product ET is equal to h. Such a reading should suggest that these objects do NOT 
admit point-like approximation since the relativity principle for free point particles requires 
straight-line uniform motion, hence, no periodicity should be allowed. 

Although the great (from pragmatic point of view) achievements of the developed theoretical 
approach, known as quantum theory, the great challenge to build an adequate description of 
individual representatives of these objects, especially of light quanta called by Lewis photons 
[19], is still to be appropriately met since the efforts made in this direction, we have to admit, 
still have not brought satisfactory results. Recall that Einstein in his late years recognizes [20] 
that "the whole fifty years of conscious brooding have not brought me nearer to the answer 
to the question "what are light quanta", and now, half a century later, theoretical physics 
still needs progress to present a satisfactory answer to the question "what is a photon". We 
consider the corresponding theoretically directed efforts as necessary and even urgent in view 
of the growing amount of definite experimental skills in manipulation with individual photons, 
in particular, in connection with the experimental advancement in the "quantum computer" 
project. The dominating modern theoretical view on microobjects is based on the notions and 
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concepts of quantum field theory (QFT) where the structure of the photon (as well as of any 
other microobject) is accounted for mainly through the so called structural function, and highly 
expensive and delicate collision experiments are planned and carried out namely in the frame of 
these concepts and methods (see the 'PHOTON' Conferences Proceedings, some recent review 
papers: [21-24]). Going not in details we just note a special feature of this QFT approach: if the 
study of a microobject leads to conclusion that it has structure, i.e. it is not point-like, then the 
corresponding constituents of this structure are considered as point-like, so the point-likeness 
stays in the theory just in a lower level. 

According to our view on PhLO we follow here an approach based on the assumption that 
the description of the available (most probably NOT arbitrary) spatial structure of photon-like 
objects can be made by continuous finite/localized functions of the three space variables. The 
difficulties met in this approach consist mainly, in our view, in finding adequate enough math- 
ematical objects and solving appropriate PDE. The lack of sufficiently reliable corresponding 
information made us look into the problem from as general as possible point of view on the 
basis of those properties of photon-like objects which may be considered as most undoubtedly 
trustful, and in some sense, identifying. The analysis made suggested that such a property seems 
to be the available and intrinsically compatible translational-rotational dynamical structure, so 
we shall focus on this property in order to sec what useful for our purpose suggestions could 
be deduced and what appropriate structures could be constructed. All these suggestions and 
structures should be the building material for a step-by-step creation of a self-consistent system. 
Prom physical point of view this should mean that the corresponding properties may combine 
to express a dynamical harmony in the inter-existence of appropriately defined subsystems of a 
finite and time stable larger physical system. 

3.2.2 The notion of photon-like object 

We begin with recalling our view that any notion of a physical object must unify two kinds 
of properties of the object considered: identifying and kinematical. The identifying properties 
being represented by quantities and relations, stay unchanged throughout the existence, i.e. 
throughout the time-evolution, of the object, they represent all the intrinsic structure and re- 
lations. The kinematical properties describe those changes, called admissible, which do NOT 
lead to destruction of the object, i.e. to the destruction of any of the identifying properties. 
Correspondingly, physics introduces two kinds of quantities and relations, identifying and kine- 
matical. From theoretical point of view the more important quantities used turn out to be the 
dynamical quantities which, as a rule, are functions of the identifying and kinematical ones, and 
the joint relations they satisfy represent the necessary interelations between them in order this 
object to survive under external influence. This view suggests to introduce the following notion 
of Photon-like object(s) (we shall use the abbreviation "PhLO" for "Photon-like object(s)"): 

PhLO are real massless time-stable physical objects with an intrinsically 
compatible translational-rotational dynamical structure. 

We give now some explanatory comments, beginning with the term real. First we emphasize 

that this term means that we consider PhLO as really existing physical objects, not as appropriate 
and helpful but imaginary (theoretical) entities. Accordingly, PhLO necessarily carry energy- 
momentum, otherwise, they could hardly be detected by physical means. Second, PhLO can 
undoubtedly be created and destroyed, so, no point-like and infinite models are reasonable: 
point-like objects are assumed to have no structure, so they can not be destroyed since there is 
no available structure to be destroyed; creation of infinite physical objects (e.g. plane waves) 
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requires infinite quantity of energy to be transformed from one kind to another during finite 
time-periods, which seems also unreasonable. Accordingly, PhLO are spatially finite and have 
to be modeled like such ones, which is the only possibility to be consistent with their "created- 
destroyed" nature. It seems hardly reasonable to believe that PhLO can not be created and 
destroyed, and that spatially infinite and indestructible physical objects may exist at all. Third, 
"spatially finite" implies that PhLO may carry only finite values of physical (conservative or 
non-conservative) quantities. In particular, the most universal physical quantity seems to be 
the energy-momentum, so the model must allow finite integral values of energy-momentum to 
be carried by the corresponding solutions. Fourth, "spatially finite" means also that PhLO 
propagate, i.e. they do not "move" like classical particles along trajectories, therefore, partial 
difi'erential equations should be used to describe their evolution in time. 

The term " massless" characterizes physically the way of propagation in terms of appropriate 
dynamical quantities: the integral 4-momentum p of a PhLO should satisfy the relation p^p^ = 0, 
meaning that its integral energy-momentum vector must be isotropic, i.e. to have zero module 
with respect to Minkowski (pseudo) metric in M^. If the object considered has spatial and 
time-stable structure, so that the translational velocity of every point where the corresponding 
field functions are different from zero must be equal to c, we have in fact null direction in the 
space-time intrinsically determined by a PhLO. Such a direction is formally defined by a null 
vector field ^, (^^ = 0. The integral trajectories of this vector field are isotropic (or null) straight 
lines as is traditionally assumed in physics, except in presence of gravity. It follows that with 
every PhLO a null straight line direction is necessarily associated, so, canonical coordinates 
"') = {x, y,z,^ = ct) on may be chosen such that in the corresponding coordinate 
frame C, to have only two non-zero components of magnitude 1: Q^^ = (0, 0, — e, 1), where e = ±1 
accounts for the two directions along the coordinate z (further such a coordinate system will be 
called ^-adapted and will be of main usage). It seems important to emphasize that our PhLO 
propagates as a whole along the (^-direction, so the corresponding energy-momentum tensor field 
Tfj,u{x, y, z, ^) of the model must satisfy the corresponding local isotropy (null) condition, namely, 
Tn^T^^'' = (summation over the repeated indices is throughout used). 

The term "translational- rotational" means that besides translational component along ^, 
the propagation necessarily demonstrates some rotational (in the general sense of this concept) 
component in such a way that both components are compatible and exist simultaneously, and 
this is an intrinsic property. It seems reasonable to expect that such kind of dynamical behavior 
should require some distinguished spatial shapes. Moreover, if the Planck relation E = hv must 
be respected throughout the evolution, the rotational component of propagation should have 
time-periodical nature with time period T = = h/E = const, and one of the two possible, 
left or right, orientations. It seems reasonable also to expect spatial periodicity of PhLO, which 
somehow to be related to the time periodicity. 

The term " dynamical structure" means that the propagation is supposed to be necessarily 
accompanied by an internal energy-m,om,entum redistribution, which may be considered in the 
model as energy-momentum exchange between (or among) some appropriately defined subsys- 
tems. It could also mean that PhLO live in a dynamical harmony with the outside world, i.e. 
any outside directed energy-momentum flow should he accompanied by a parallel inside directed 
energy-momentum flow. 

Finally, note that if the time periodicity and the spatial periodicity should be consistent 
with each other somehow, the simplest integral feature of such compatability would seem like 
this: the spatial size A along the translational component of propagation is equal to cT: A = cT, 
where A is some finite positive characteristic constant of the corresponding solution. This would 
mean that every individual PhLO determines its own length/time scale. 
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It is important to note now the following. We don't know what mathematical objects are 
appropriate for describing PhLO, so, our first task is to come to such mathematical objects 
having in view what we mean under PhLO. The next Section is devoted namely to find math- 
ematical structures that are adequate enough to the above introduced notion for PhLO and 
earring rich enough fiexability to meet all requirements for a field theory of spatially finite and 
time-stable physical objects with dynamical structure. Our hope is that the ideas and concepts 
connected with the Frobenius integrability theory seem to represent the most adequate part of 
mathematics for this purpose. 

3.3 Curvature of Distributions and Physical Interaction 

3.3.1 The general idea for geometrization of local physical interaction 

We begin with a short motivation for this choice of mathematics directed to the readers already 
acquanted with Frobenius integrability theory, and right after this we shall carefully introduce 
the necessary mathematics. 

Any physical system with a dynamical structure is characterized by some internal energy- 
momentum redistributions, i.e. internal energy-momentum fluxes, during evolution. Any time- 
stable compatible system of energy- momentum fluxes (as well as fluxes of other interesting for the 
case physical quantities subject to change during evolution, but we limit ourselves just to energy- 
momentum fluxes here) can be considered mathematically as generated by a compatible system 
of vector fields. A physically isolated and interelated time-stable system of energy-momentum 
fluxes can be considered to correspond directly or indirectly to a completely integrable distribu- 
tion A of vector fields (or differential system [25]) according to the principle: some local objects 
can generate integral object. Every nonintegrable distribution on a manifold defines its own cur- 
vature form (given further in the section). Let Ai and A2 be two nonintegrable distributions on 
the same manifold with corresponding curvature forms Oi and ^^2, each of them carries couples 
of vector fields inside their distributions outside Ai and A2 correspondingly, i.e. Oi(Yi,y2) / 
is out of Ai and $12(^1,^2) 7^ is out of A2, where {Yi,Y2) live in Ai and (^1,^2) live in 
A2. Let now Ai and A2 characterize two locally interacting physical systems, or two locally 
interacting subsystems of a larger physical system. It seems reasonable to assume as a workong 
tool the following geometrization of the concept of local physical interaction: two nonintegrable 
distributions Ai and A2 on a manifold will be said to interact infinitesimally (or locally) if some 
of the nonzero values of the corresponding two curvature forms ^1/0,2 live respectively in A2/A1. 

The above geometric concept of infinitesimal interaction is motivated by the fact that, in gen- 
eral, an integrable distribution A may contain various nonintegrable subdistributions Ai, A2, . . . 
which subdistributions may be associated physically with interacting subsytems of a larger time 
stable physical system. Any physical interaction between 2 subsystems is necessarily accom- 
panied with available energy-momentum exchange between them, this could be understood 
mathematically as nonintegrability of each of the two subdistributions of A and could be nat- 
urally measured directly or indirectly by the corresponding curvatures. For example, if A is 
an integrable 3-dimensional distribution represented by the vector fields {Xi,X2,X-i) then we 
may have, in general, three non-integrable, i.e. geometrically interacting, 2-dimensional subdis- 
tributions (Xi,X2), {Xi,X^), (X2, Xa). Finally, some interaction with the outside world can be 
described by curvatures of distributions (and their subdistributions) in which elements of A and 
vector fields outside A are involved (such processes will not be considered in this paper). 

The above considerations launch the general idea to consider the concept of Frobenius cur- 
vature as a natural and universal mathematical tool for describing local physical interaction be- 
tween/among the relatively stable subsystems of the physical world. In other words, the Frobe- 
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nius curvature appears as appropriate mathematical tool describing formally the 
possible ability two continuous systems to recognize each other as physically inter- 
acting partners. 

Two formal aspects of the obove idea exist. The first applies directly the Frobenius integra- 
bility machinary [25] , while the second one (been developed recently) is known as nonlinear 
connections [26]. We consider now briefly the first one. 

3.3.2 Frobenius integrability, curvature and local physical interaction 

A p-dimensional distribution Ap on a n-dimensional manifold M" is defined by associating to 
each point x G M" a p-dimensional subspace of the tangent space at this point: c T^-M", x G 
M",l < p < n. Let the system of vector fields {Xi,X2, . . . ,Xp} represent this distribution, 
so {Xi{x),X2{x),. . . ,Xp{x)}, X G M"^, I < p < n, satisfy Xi{x) A X2{x) A ...,AXp{x) ^ 
0, X E M", and represent a basis of A§. According to the Frobenius integrability theorem 
(further all manifolds are assumed smooth and finite dimensional and all objects defined on 
M" are also assumed smooth) Ap is completely integrable, i.e. through every point x G M" 
passes a p-dimensional submanifold such that all elements of Ap are tangent to N^, iff 
all Lie brackets , i, j = 1,2, . . . ,p, are representable linearly through the very Xi,i = 

1,2,... ,p : [Xi,Xj] = C^jXk, where C^j are functions. Clearly, an easy way to find out if a 
distribution is completely integrable is to check if the exterior products 

[Xi,Xj]AXi{x)AX2{x)A...,AXp{x),xeM''; i, j = 1,2, . . . ,p (3.1) 

are identically zero. If this is not the case (which means that at least one such Lie bracket 
"sticks out" of the distribution Ap) then the corresponding coefficients, which are multilinear 
combinations of the components of the vector fields and their derivatives, represent the corre- 
sponding curvatures. We note finally that if two subdistributions contain at least one common 
vector field it seems naturally to expect interaction. 

In the dual formulation of Frobenius theorem in terms of differential 1-forms (i.e. Pfaff 
forms), having the distribution Ap , we look for (n — p)-Pfaff forms {a^,a^, . . . ,a""~^), i.e. a 
(n — p)-codistribution A* , such that {a"^,Xj) = 0, and a^{x) Aa^{x) A ■■■ A a""P(x) / 0, 
m = 1,2, . . . ,n — p, j = 1,2, . . . ,p,x e M". Then the integrability of the distribution Ap is 
equivalent to the requirements 

da"^ Aa^ Aa^ A--- Aa"~P = 0, m = 1,2, . . . , {n - p), (3.2) 

where d is the exterior derivative. 

Since the idea of curvature associated with, for example, an arbitrary 2-dimensional distribu- 
tion {X, Y) is to find out if the Lie bracket [X, Y] (p) has components along vectors outside the 
2-plane defined by (Xp,Yp), in our case we have to evaluate the quantities (a™, [X, y]), where 
all linearly independent 1-forms annihilate {X,Y) : (q;"^,X) = {a"^,Y) = 0. In view of the 
formula 

da^{X,Y) = X{{a^,Y)) - Y{{a^,X)) - {a"',[X,Y]) = -{a"',[X,Y]) 

we may introduce explicitly the curvature 2-form for the distribution A(X) = (Xi, . . . , Xp). 
In fact, if A(y) = (Yi, . . . ,Yn-p) define a distribution which is complimentary (in the sense of 
direct sum) to A{X) and = 0, = (5^, i.e. {Yi, . . . ,F„-p) and (a\ . . . ,a"-P) 

are dual bases, then the corresponding curvature 2-form J1a(x) should be defined by 

^A{X) = -dcT ® Ym, since Qa(x) {Xi,Xj) = -da"*(Xi, Xj)Ym = (a"^, [X^, Xj])Ym, (3.3) 
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where it is meant here that ^a(x) is restricted to the distribution {X^, . . . ,Xp). Hence, if we 
call the distribution {Xi, . . . , Xp) horizontal and the complimentary distribution (Yi, . . . , Yn-p) 
vertical, then the corresponding curvature 2-form acquires the status of vertical bundle valued 
2-form. We see that the curvature 2-form distinguishes those couples of vector fields inside A(X) 
the Lie brackets of which define outside A(X) directed flows, and so, do not allowing to find 
integral manifold of A{X). Clearly, the supposition here for dimensional complementarity of the 
two distributions A(X) and A(Y) is not essential for the idea of geometrical interaction, i.e. the 
distribution A(F) ^ A(X) may be any other distribution on the same manifold with dimension 
smaller than {n — p), so that m = 1,2,. . . ,q < {n — p) in general, the important moment is 
that the two distributions (or subdistributions) can " communicate" differentially through their 
curvature 2-forms. 

Hence, from physical point of view, if the quantities Xj) are meant to be used for 

building the components of the energy- momentum locally transferred from the system A(X) to 
the system A{Y), then, naturally, we have to make use of the quantities ^A{Y)0^myYn) to build 
the components of the energy-momentum transferred from A(Y) to A(X). 

It deserves to note that this formalism allows a dynamical equilibrium between the two 
systems A(y) and A{X) to be described: each system to gain from the other as much energy- 
momentum as it loses, and this to take place at every space-time point. Therefore, if W(^x,Y) 
denotes the energy-momentum transferred locally from A(X) to A(y), W(^y,x) denotes the 
energy-momentum transferred locally from A{Y) to A{X), and SW(^x) a^nd 6W(y) denote respec- 
tively the local energy-momentum changes of the two systems A{X) and A{Y), then according 
to the local energy-momentum conservation law we can write 

SW(^X) = W(^Y,X) + W(^x,Y)^ ^^(Y) = -(.W(^X,Y) + W{Y,X)) = -SW(^X), 

i.e. A{X) and A(Y) are phisically compatible, or able to interact. For the case of dynamical 
equilibrium we have W(^x,Y) ~ ~'^{Y,x): so in such a case we obtain 

5W^x)=0, SW(^y) = 0, W(^Y,x) + W^x,Y)=0- (3.4) 

As for how to build explicitly the corresponding representatives of the energy-momentum fluxes, 
probably, universal procedure can not be offered. The most simple procedure seems to be 
to "project" the curvature values ilA{x)i^ij ^j) ^A{Y)0^m,Yn) on the corresponding co- 
distribution volume forms, i.e. to consider the corresponding inner products i{^{Xi, Xj)){a^ A 
A • • • A a"'~P). For every special case, however, appropriate quantities constructed out of the 
members of the introduced distributions and co-distributions must be worked out. 



3.4 PhLO Dynamical Structure in Terms of Frobenius 
Curvature 

We consider the Minkowski space-time M = (M"^,?]) with signature sign{r]) = (— ,—,—,-!-) 

related to the standard global coordinates (x^, x^, x^, x"^) = {x,y,z,^ = ct),the natural volume 
form cOo = ^/\r]\dx^ A dx^ A dx^ A dx* = dx A dy A dz A d^, and the Hodge star * defined by 
a A*P = —r]{a, (3)uJo. 

In view of our concept of PhLO we introduce the null vector field C, (^^ = 0, which in the 
C-adapted coordinates (throughout used further) is assumed to look as follows: 
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Let's denote the corresponding to ( completely integrable 3-dimensional Pfaff system by A*(i^). 
Thus, A*{() can be generated by any three linearly independent 1-forms (ai, 0:2, 0:3) which 
annihilate (, i.e. 

Q!i(C) = "2(0 = «3(C) = 0; Oil A a2 A a3 0. 
Instead of (ai, 02, 03) we introduce the notation {A, C) and define ( to be the 77-corresponding 
1-form to ^: 

C = edz + d^, so, (C,C)=0, (3-6) 
where (, ) is the coupling between forms and vectors. 

Now, since is closed, it defines 1-dimensional completely integrable Pfaff system, so, we 
have the corresponding completely integrable distribution {A,A*,() : {C,A) = {C,A*) = 0. We 
shall restrict our further study to PhLO of electromagnetic nature according to the following 

Definition: We shall call a PhLO electromagnetic if the following conditions hold: 

1. the vector fields (A, A*) have no components along ^, 

2. (A, A*) are ry-corrcsponding to (^4, ^*) respectively . 

3. {A, A*) = 0, {A, A) = {A*, A*) . 

Remark. These relations formalize knowledge from Classical electrodynamics (CED). In fact, 
our vector fields {A, A*) are meant to represent what we call in CED electric and magnetic 
components of a free time-dependent electromagnetic field, where, as we have mentioned several 
times, the translational propagation of the field energy-momentum along a fixed null direction 
with the velocity "c" is possible only if the two invariants Ii = — and I2 = 2E.B are 
zero, because only in such a case the electromagnetic energy-momentum tensor T^jy satisfies 
T^^T^'^ = and has unique null eigen direction. So it seems naturally to consider this property 
as intrinsic for the field and to choose it as a starting point. Moreover, in such a case the relation 
(Ii)^ -I- (12)^ = is equivalent to E-^ -|- = 2|E x B| and this relation shows that this is the 
only case when the field momentum can not be made equal to zero by means of frame change. 
Together with the fact that the spatial direction of translational encrgy-momcntTim propagation 
is determined by E x B, this motivates to introduce the vector field Q in this form and to assume 
the properties 1-3 in the above definition. 

Prom the above conditions it follows that in the ^-adapted coordinate system we have 

A 1 , A* 1 1 A d d d d 

A = udx + pdy, A = —epdx + eudy; A = —u— Pt;-, A* = £p— en—, 

ox ay ox Oy 

where e = ±1, and (u,p) are two smooth functions on M. 

The completely integrable 3-dimensional Pfaff system {A, A* , Q contains three 2-dimensional 
subsystems: {A, A*), {A,() and {A*,(). We have the following 

Proposition 1. The following relations hold: 

dAAAAA* = 0; dA* aA* aA = 0; 

dAAAA( = e[w(pg - epz) - p{u^ - en^)] Wo; 
dA* A A* AC, = £ \u{p^ — epz) — p{u^ — £Uz)]oOo- 
Proof. Immediately verified. 

These relations say that the 2-dimensional Pfaff system {A, A*) is completely integrable for any 
choice of the two functions (u,p), while the two 2-dimensional Pfaff systems {A,^) and (^*,C) 
are NOT completely integrable in general, and the same curvature factor 

R = u{p^ - epz) - p{u^ - euz) 
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determines their nonintegrability. 

Correspondingly, the 3-dimensional completely integrable distribution (or differential system) 
A((^) contains three 2-dimensional subsystems: {A, A*), (^,C) ^-nd iA*,C)- We have the 

Proposition 2. The following relations hold (recall that [-'^j i^] denotes the Lie bracket): 

[A, A*]aAaA* = 0, (3.7) 

d d 
[A, C] = {u^ - eu,)— + {jp^ - ep,) — , (3.8) 

[A*X] = -e{p^-ep,)—+£{u^-£u,) — . (3.9) 
Proof. Immediately verified. 

From these last relations (3.7-3.9) it follows that the distribution {A^A*) is completely in- 
tegrable, and it can be easily shown that the two distributions (A, ^) and (A*,^) would be 
completely integrable only if the same curvature factor 

R = u(jp^ - epz) - p{u^ - £Uz) (3.10) 

is zero (the elementary proof is omitted). 
As it should be, the two projections 

{A, [A*,C]) = -{A*, [A, (]) = eu{p^ - ep,) - ep{u^ - eu,) = -eR 

are nonzero and give (up to a sign) the same factor R. The same curvature factor appears, of 
course, as coefficient in the exterior products [A* , (^] A A* A C and [^4, A A A C. In fact, we obtain 

[A-*,C]AA-*AC = -[AC-]AAAC = -eR^A|A|+RAA|A|. 

On the other hand, for the other two projections we obtain 

(AKCl) = (^*,[^*,C]) = \[{u'+p\-e{v?+p\]. (3.11) 
Clearly, the last relation (3.11) may be put in terms of the Lie derivative as 

^L-^iu-" +p^) = -h^{A,A) = -{A,L^A) = -{A*,L^A*). 

Remark. Further in the paper we shall denote y^ti^ + p'^ = (p. 
Proposition 3. There is a function ipiu^p) such, that 



Proof. It is immediately verified that ^ = arctan ^ is such one. 

We note that the function ip has a natural interpretation of phase because of the easily 
verified now relations u = 0cos'(/', P = (psinip, and (/) acquires the status of amplitude, i.e. 
energy density. Since the transformation {u,p) — > (^, V') is non-degenerate this allows to work 
with the two functions ((?!>, V') instead of {u,p). 
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Prom Prop.3 we have 



R = ,^L,i, = - s,.) L,^ = ^ = (3.12) 



where T(d^, d^) is the coordinate-free definition of the energy density. 

This last formula (3.12) shows something very important: at any 7^ the curvature R 
will NOT be zero only if L^ip 7^ 0, which admits in principle availability of rotation. In fact, 
lack of rotation would mean that (j) and arc running waves along ^. The relation L^ijj 7^ 
means, however, that rotational properties are possible in general, and some of these properties 
are carried by the phase ^. It follows that in such a case the translational component of 
propagation along ( (which is supposed to be available) must be determined essentially, and 
most probably entirely, by (p. In particular, we could expect the relation L^<p = to hold, and 
if this happens, then the rotational component of propagation will be represented entirely by 
the phase V", and, more specially, by the curvature factor R 7^ 0, so, the objects we are going to 
describe may have compatible translational-rotational dynamical structure. Finally, (3.12) may 
be considered as a definition for the phase function tp. 

We are going now to represent some relations, analogical to the energy-momentum relations 
in classical electrodynamics, determined by some 2-form F, in terms of the Probenius curvatures 
given above. 

The two nonintegrable Pfaff systems (^4, Q and {A*, () carry two volume 2- forms: 

G = AaC and G* =A* AC, 

and the two corresponding distributions define the 2-vectors 

G = AaC, and G* = A* A (. 

Making use now of the Hodge *-operator, we can verify the relation: G* = *G. Now G and G* 
define the (l,l)-tensor, called stress-energy-momentum tensor T^, according to the rule 

T;: = -^[G,,G'''^ + {G*u{G*r], 

and the divergence of this tensor field can be represented in the form 

V.r;=[z(G')dG]^+[i(G*)dG*]^, 

where G and G* coincide with the metric-corresponding contravarint tensor fields, and i{G) = 
i{Q o i{A), i{G*) = i{Q o i{A*), i{X) is the standard insertion operator in the exterior algebra 
of differential forms on defined by the vector field X. So, we shall need the quantities 

i{G)dG, i{G*)dG*, i{G*)dG, i{G)dG*. 

Having in view the explicit expressions for A, A*,(, A, A* and ( we obtain 

i{G)dG = i{G*)dG* = ^L^ (</.2) . C , (3.13) 

also, we obtain 

i{G*)dG = -i{G)dG* = 

u{p^ - epz) - p{u^ - euz) dz + e u{p^ - ep^) - p{u^ - eu^) = eRC. (3.14) 
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If F and H are correspondingly 2 and 3 forms on M we have the relation 



*{F A *H) = i{F)H = F''" Hi^^^dx" , n<v. 

Therefore, since G* = *G, 

z(G*)dG = -*((5*GA*G), -i{G)dG* = -*{5G^G), so, S*G A*G = SG AG = eK*C. (3.15) 

In the following formulae we must keep in mind the relations d(" = 0, {A, A*) = {A*, A) = 
{(^A*) = iCA) = 0, and {Af = {A*f = {A, A) = {A*, A*) = -{u^+p"^) = -^^ = _|^|2 = 

-|^*|2 = -|i|2 = -|^*|2. 

In view of these formulae and the required duality in the definition of the curvature form 
(3.3), the two distributions (^, C) aiid (^*,C) determine the following two curvature forms ft 
and 9,*: 

^ -A* A* ^ A* A* -A A _^A A 

\A*\ \A*\ \A*\ \A*\ \A\ \A\ \A\ \A\ 

Denoting Zq, = n{A,C), = n{A*,C), Zq* = n*{A,C) and Z^, = n*{A*,Q we obtain 

c-T? _ A* A pTt 

Zn = -^A*, Z*^ = -—^L-^{<^% Z^. = -—^L-^{<^% = 0^^" (^'l^) 

The following relations express the connection between the curvatures and the energy-momentum 
characteristics. 

z(Zn)(AAC) = 0, i{Za){A* AC,)=eIi.C,=-i{G)dG* = i{G*)dG, (3.17) 
i{ZQ*){A* aQ = {), i{Z'^.){AACi = -eK.(: = i{G)dG* = -i{G*)dG, (3.18) 



i{Z-^){A A C) = 0, i{Z*^){A* A C) = \Lz{ct>'').C = i{G)dG = i{G*)dG\ (3.19) 

1 



2 



i{Z*^,){A* A C) = 0, i{Z^*){A A C) = ^L-^m-C = i{G)dG = i{G*)dG\ (3.20) 
It follows from these relations that in case of dynamical equilibrium we shall have 
L^-(02) = 0, i{G)dG = 0, i{G*)dG* = 0, i{G*)dG + i{G)dG* = 0. 

Resuming, we can say that Probenius intcgrability viewpoint suggests to make use of one 
completely integrable 3-dimensional distribution (rcsp. Pfaff system) consisting of one isotropic 
and two space-like vector fields (resp. 1-forms), such that the corresponding 2-dimensional spa- 
tial subdistribution {A, A*) (resp. Pfaff system {A, A*)) defines a completely integrable system, 
and the rest two 2-dimensional subdistributions {A,C,) and (^*,C) (resp. Pfaff systems (^, C) 
and {A*,C,)) are NON-integrable in general and give the same curvature. This curvature may 
be used to build quantities, physically interpreted as energy-momentum internal exchanges be- 
tween the corresponding two subsystems {A,Q and (^*,C) (resp.(^, (^) and {A*,C,)). Moreover, 
rotational component of propagation will be available only if the curvature R is nonzero, i.e. 
only if an internal energy-momentum exchange takes place. Wc sec that all physically important 
characteristics and relations, describing the translational and rotational components of propa- 
gation, can be expressed in terms of the corresponding Probenius curvature. We'll see that this 
holds also for some integral characteristics of PhLO. 
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3.5 PhLO Dynamical Structure in Terms of Non- linear 
Connections 



3.5.1 Projections and algebraic curvatures 



The projections are linear maps P in a linear space W"" (under linear space we mean here 
module over a ring, or vector space over a field) sending all elements of to some subspace 
P(W^") C W, such that P o P = P. Let (ei, . . . , e^, . . . , en) and {e^, . . . ,£P, . . . be two 
dual bases: < e^, >= 5^, ^u, = 1, . . . , n, and let Nf, i = 1, . . . ,p ; a = p + 1, . . . ,n he the 
corresponding to P [p x (n — p)] matrix of rank (n — p). We define another couple of dual bases: 



k^ = {ei + Ntea,ea); co" = {e\ a' - N^e^), j 



,p; a,b = p+l,...,n. 



Now the identity map idw^ 
idwn = Lo^ k,y = Lo^ 



®ky acquires the form 



ki + oo ®kb = 



{ei + Ntea) + {e'-N'^e^] 



(3.21) 



Wc obtain two projections: Py 



Nh^) 66 and Ph 



(cj + Nfea) such that 



KerPv = ImPn and KerPn = ImPy, also, Ph = idwn — Py- Hence, IF" = KerPy (B ImPy = 
KerPH^BlmPn- Usually Py is called vertical projection, and Ph is called horizontal projection. 

Let now (f) and i/j be two arbitrary linear maps in a module 971, 58:SDTx9Jt— >^9Jtbea 
binar map satisfying *B(x + z, y) = QS(x, y) + QS(z, y) and 5S(x, y + z) = QS(x, y) + *8(x, z), and 
(x, y, z) be three arbitrary elements of 93t. We consider the expression 



^(!B;0,V)(x,y) 



!B(</>(x), V(y)) + 5B(V^(x), 0(y)) + </) o V;(5B(x, y)) + V o <^(iB(x, y)) 



-0(©(x,^(y))) - </.(©(V'(x),y)) - V(»(x,</.(y))) - V(»(<^(x),y)) . 
Assuming (j) = ip arc projections in dJl denoted by P, this expression becomes 

^(<8;P)(x,y) = P(«8(x,y))+«8(P(x),P(y)) -P(«(x,P(y))) -P(»(P(x),y)) . 

Denoting the identity map of 9Jl by id and adding and subtracting P *B(P(x), P(y)) , after 
some elementary transformations we obtain 

Ai^;P) (x, y) = p[3 [{id - P) (x), {id - P) (y)] ]+{id- P) [P(x) , P(y] 

Recalling that P and {id — P) project on two subspaces of the direct sum of which generates 
9Jt, and naming P as vertical projection denoted by V, then {id — P), denoted by H, gets 
naturally the name horizontal projection. So the above expression gets the final form of 



^(58;P)(x,y) = Ffe[/f(x),i7(y)] 



+ H 



« [F(x), y(y)] J = 7^p(©; X, y) + 1lp{^- x, y). 

(3.22) 

Hence, the first term on the right, 7^p(Q5;x,y), which may be called '^-algebraic curvature of P, 
measures the vertical component of the ©-image of the horizontal projections of (x, y), and then 
the second term '^p(®i;x, y), acquiring the name of ^-algebraic cocurvature of P, measures the 
horizontal component of the IB-image of the vertical projections of (x, y). 

We carry now this pure algebraic construction to the tangent bundle of a smooth manifold 
M", where the above binar map 53 will be interpreted as the Lie bracket of vector fields, 
and the linear maps will be just linear endomorphisms of the tangent/cotangent bundles of 
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M". Under these assumptions the quantity A{^,^) is called Nijenhuis bracket of the two 
linear endomorphisms $ and ^, and is usually denoted by [^*, ^']. It has two important for us 
properties: the first one is that ^] is linear with respect to the smooth functions on the 
manifold, so, the Nijenhuis bracket allows, starting with two (1, l)-tensors on M'", to construct 
through differentiations a 2-form that is valued in the tangent bundle of M'*; the second property 
is that if $ = * then is not necessarily zero. 

3.5.2 Nonlinecir connections 

Let now {x^,...,x'^) be any local coordinate system on our real manifold M". We have the 
corresponding local frames {dx^ , . . . , dx"^} and {d.j.i , . . . , dx"} ■ Let for each x G M we are 
given a projection of the same constant rank (n — p), i.e. p does not depend on x, in every 
tangent space Tx{M). The space Ker{Px) C Tx{M) is usually called P-horizontal, and the space 
Im{Px) C Tx{M) then is called P-vertical. Thus, we have two distributions on M the direct sum 
of which gives the tangent bundle: T{M) = Ker{P)® Im{P). The above algebraic construction 
shows that each of these two distributions can be endowed with corresponding 2-form, valued in 
the other distribution, and depending on the same binar operation in TM". As we mentioned, 
the choice Si= Lie bracket leads to tensor field. Therefore, assuming this choice, we say that P 
defines a nonlinear connection on M. Denoting by TZ the so defined curvature 2-form of P and 
by TZ the corresponding cocurvature 2-form of P, by Vp and Hp the corresponding vertical and 
horizontal prjections, we can write 

[P,P]{X,Y) = n{X,Y) + n{X,Y), (3.23) 

where 

7^(X, Y) = Vp{[HpX, HpY]) , iiiX, Y) = Hp{[VpX, VpY]) , 

(X, Y) are any two vector fields and the Lie bracket is denoted by [, ] . Recalling the contents of 
the preceding section, it can be shown that TZ{X, Y) ^ measures the nonintcgrability of the 
corresponding horizontal distribution, and TZ(X, Y) ^ measures the nonintcgrability of the 
corresponding vertical distribution. 

If the vertical distribution is given before-hand and is completely integrablc, i.e. TZ = 0, then 
TZ{X, Y) is called curvature of the nonlinear connection P if there exist at least one couple of 
horizontal vector fields (X, Y) such that TZ{X, Y) ^ 0. 

3.5.3 Photon-like nonlinecir connections 

We assume now that our manifold is endowed with standard coordinates [x^ 
x,y,z,(^ = ct), and make some preliminary considerations in order to make the choice of our 
projection P : TW^ TW^ consistent with the introduced concept of PhLO. The intrinsically 
defined straight-line translational component of propagation of the PhLO will be assumed to 
be parallel to the coordinate plane (z,^). Also, ^ and ^ will be vertical coordinate fields, 

so every vertical vector field Y can be represented by Y — u P '§yi where (u, p) are two 
functions on M^. It is easy to check that any two such linearly independent vertical vector 
fields Yi and Y2 define an integrable distribution, hence, the corresponding curvature will be 
zero. It seems very natural to choose Yi and I2 to coincide correspondingly with the vertical 
projections P{-§^) and -P(^)- Moreover, let's restrict ourselves to PhLO of electromagnetic 
nature and denote further the verical projection by V . Then, since this vertical structure is 
meant to be smoothly straight-line translated along the plane (-2,^) with the velocity of light. 
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a natural suggestion comes to mind these two projections Yi = V{^) and Y2 = V{^) to be 
physically interpreted as representatives of the electric and magnetic components. Now we know 
from classical electrodynamics that the situation described corresponds to zero invariants of the 
electromagnetic field, therefore, we may assume that Yi and Y2 are ortogonal to each other and 
with the same modules with respect to the euclidean metric in the 2-dimensional space spent by 
^ and It follows that the essential components of Yi and Y2 should be expressible only with 
two independent functions {u,p). The conclusion is that our projection should depend only on 
{u,p). Finally, we note that these assumptions lead to the horizontal nature of dz and d^. 

Note that if the translational component of propagation is along the vector field ^ then we 
can define two new distributions : {Yi,() and {Y2,C), which do not seem to be integrable in 
general even if ( has constant components as it will be in our case. Since these two distributions 
are nontrivially intersected (they have a common member (), it is natural to consider them 
as geometrical images of two interacting physical subsystems of our PhLO. Hence, we must 
introduce two projections with the same image space but with different kernal spaces, and the 
components of both projections must depend only on the two functions {u,p). 

Let now {u,p) be two smooth functions on and e = ±1 . We introduce two projections 
V and V in TM^ as follows: 



V = dx(Si-^ + dy®-^ eudz® ^ udi® ^ epdz® ^ pdf®^, (3.24) 



d 

dx 


+ dy^ 


d 
dy 


— £udz ( 


9 

^ — -udC(^ 
dx 


9 

^ epdz 

dx 


9 

® pd^^ 

dy 


d 
dy 


d 
dx 


+ dy<i 


d 
dy 


+ pdziSi 


■^ + epd^^ 
dx 


i — udz<Si 

dx 


eud^^ 

dy 


d 
dy 



(3.25) 

ux ux uy uy 

So, in both cases we consider (^, ^) as vertical vector fields, and {dz, d^) as horizontal 1-forms. 

By corresponding transpositions we can determine projections V* and V* in the cotangent 
bundle r*M4. 

1 9 1 d ,9 , d , d ,9 

V = dx® — +dy ® eudx® udx ® — - epdy ® pdy®—, 

dx dy dz dt, dz d$, 

■Tr* , 9 ^ d ,9 ,9 ,9 ,9 

V =dx® — + dy® — + pdx® — + epdx® — - udy<i^- eudy ® —. 

dx dy dz d^ dz dt, 

The corresponding horizontal projections, denoted by {H,H;H*H*) look as follows: 

d d d d d d 

H = dz ® — + d^,® — + eudz ® — + ud^® — + epdz ® — +pd^ ® — , 
dz dt, dx dx dy dy 

~ d d d d d d 

H = dz® — + d^® — - pdz®- epd^® — + udz® — + eud^® —, 

dz d^ dx dx dy dy 

TT* 1 9 d , d ,9 , d ,9 

H = dz ® — + dt,® — + eudx ® — + udx ® — + epdy ® — + pdy ® 

dz d^ dz d^ dz d^ 

^T* , 9 ,^ d , d ,9 ,9 ,9 

H = dz® — +d^® — -pdx ® epdx ® + udy® — + eudy ® —. 

dz dt, dz d^ dz dt, 

The corresponding matrices look like: 



V 



1 





—su 


—u 





1 


—ep 


-p 



























H 









eu 


u 








ep 


p 








1 














1 
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V* = 



1 




























1 
































, H* = 











—£U 


—ep 




eu 


ep 


1 


—u 


-P 










u 


P 





1 



v = 



1 





p 


ep 





1 


—u 


—eu 



























H = 









-p 


—ep 








u 


e u 








1 














1 



V* = 



1 




























1 
































, H* = 











p 


—u 


-p 


u 


1 


ep 


—eu 










—ep 


eu 
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The projections of the coordinate bases are: 

\dx' dy^ dz' ) ' \dx^ dy^ ^^dx ^'^dy' ^ dx ^dy)^ 



d d d d d d 
dx'dy'dz'd^J V ' ' ^ dx ^ dy dz^^dx ^ dy 

{dx, dy, dz, d$,) .V* = {dx — eudz — ud^, dy — epdz — pd^, 0, 0) 
{dx, dy, dz, d^) .H* = [eudz + ud^, epdz + pd^, dz, d^) 

dddd\^fddd d d d 

dx^ dy^ dz^ d^j \dx^ dy'^ dx dy^ ^ dx dy 

d d d d\ ~ f d d d d d d 

dx dy dz dt,) \ dx dy dz dx dy d^ 

{dx, dy, dz, d^) .V* = {dx + pdz + epd^, dy — udz — eud^^, 0, 0) 

{dz, d^, dx, dy) .H* = {—pdz — epd^,udz + eud^, dz, d$,) . 

We compute now the two curvature 2-forms TZ and TZ. The components 7?.^^ of TZ in coordinate 

basis are given by Vp (j^H , H -^Y^ , and the only nonzero components are just 

= ^34 = -^Uz), Tll^ = 7^i4 = -e{p^ - epz). 

For the nonzero components of TZ we obtain 

= ^34 = {Pi - spz), = nli = -{h-^ ^^)- 

The corresponding two curvature forms are: 

TZ = —e{u^ — euz)dz /\d£,® ^ e{jp(^ — epz)dz A (g) ^ (3.26) 

TZ = {p^- epz)dz Ad^^ — -{u^- euz)dz A (g) — . (3.27) 
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We 



V 



obtain (in our coordinate system): —\tr {V o H*) = —^tr (v o H*^ = +p^, and 



H 



d_ 

dz 



d_ 



H 



d_ 

dz 



M 



d_ 



V 



dz 



dz 



-£{u^ -£u,)^- e{p^ - £Pz)^ = Zi, 
^ (P« - ep,)-^ - {u^ - su,)^ = Z2, 



where Zi and Z2 coincide with the values of the two curvature forms TZ and TZ on the coordinate 
vector fields ^ and ^ respectively: 



Zi =7^ 



d_ d_ 



Zo = n 



d_ d_ 



We evaluate now the vertical 2-form V*{dx) A V*{dy) on the bivector Zi A Z2 and obtain e/C^, 
where 



{u^ -euz) + {p^ - epz 



An important parameter, having dimension of length (the coordinates are assumed to have 
dimension of length) and denoted by Iq-, turns out to be the square root of the quantity 



-\tr{VoH* 



V? + p^ 



+ ip^- epz 



^2• 



Clearly, if Iq is finite constant it could be interpreted as some parameter of extension of the 
PhLO described, so it could be used as identification parameter in the dynamical equations and 
in lagrangians, but only if {u^ — euz) 7^ and {p^ — epz) / 0. This goes along with our concept 
of PhLO which does not admit spatially infinite extensions. Finally we'd like to note that the 
right-hand side of the above relation does not depend on which projection or y is used, i.e. 
[V*{dx) A V*{dy)]{Zi A Z2) =eK? too, so 



1^ 



triVoH* 



-\tr[V oH* 



/C2 



(3.28) 



The parameter Iq has the following symmetry. Denote by = dx (^i) + dy <S) then 

V = Vo + Vi and V = Vo + Vi, where, in our coordinates, Vi and Vi can be seen above how they 
look like. We form now W = aVi — hVi and W = hV\ + aV\, where (a, h) are two arbitrary real 
numbers. The components of the corresponding linear maps Pw = Vo + W and = Vo + W 

can be obtained through the substitutions: u {au + ehp); p [ehp — ap), and, obviously, Pw 
and P^y are projections. Now, —^tr{VoH*) transforms to (a^ + h'^){u^ +P^) and K? transforms 
to (a^ + 6^)[(m^ — euzf' + (p^ — ep^)^], so, lo{V, V) = lo(W, W). This corresponds in some sense to 
the dual symmetry of classical vacuum electrodynamics. We note also that the squared modules 
of the two curvature forms |7?.p and \'R,\^ are equal to {u^ — euz)"^ + {p^^ — £Pz)^ in our coordinates, 
therefore, the nonzero values of \R.\^ and as well as the finite value of Iq guarantee that 

the two functions u and p are NOT plane waves. Finally, the phase function may be defined by 
the relations 

. , <A,Z2> <A*,Zi> 



-\tr{VoH*) 
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3.5.4 Electromagnetic PhLO in terms of non-lineeir connections 

Recall that the relativistic formulation of classical electrodynamics in vacuum [p = 0) is based 
on the following assumptions. The configuration space is the Minkowski space-time M = (M^, 77) 
where r] is the pseudomctric with sign{r]) = (— , — , — , +) with the corresponding volume 4-form 
LJo = dx A dy A dz A d^ and Hodge star * defined by a A */3 = —rj{a, j3)u)o- The electromagnetic 
filed is describe by two closed 2-forms (F, *F) : dF = 0, d * F = 0. The physical characteristics 
of the field are deduced from the following stress-energy-momentum tensor field 

r/(F, *F) = [F^.F-'^ + i*F)^,i*Fr] . (3.29) 

In the non-vacuum case the allowed energy-momentum exchange with other physical systems is 
given in general by the divergence 



V T" — - 



F"^(dF),^^ + {*Fr^{d * F)^^^ = F^,{5Fy + (*F)^,(5 * F)^ (3.30) 



where 5 = *d* is the coderivative. If the field is free: dF = 0, d * F = 0, this divergence is 
obviously equal to zero on the vacuum solutions since its both terms arc zero. Therefore, energy- 
momentum exchange between the two partner-fields F and *F, which should be expressed by the 
terms (*F)°^^(dF)«;3^ and F"^(d * F)^^^ is NOT allowed on the solutions of dF = 0, d * F = 
0. This shows that the widely used 4-potential approach (even if two 4-potentials A, A* are 
introduced so that dA = F, dA* = *F locally) to these equations excludes any possibility 
to individualize two energy-momentum exchanging time-stable subsystems of the field that are 
mathematically represented by F and *F. 

On the contrary, as we have mentioned several times, our concept of PhLO does NOT 
exclude such two physically interacting subsystems of the field to really exist, and therefore, 
to be mathematically individualized. The intrinsically connected two projections V and V 
and the corresponding two curvature forms give the mathematical realization of this idea: V 
and V individualize the two subsystems, and the corresponding two curvature 2-forms TZ and 
TZ represent the instruments by means of which the available mutual local energy-momentum 
exchange between these two subsystems could be described. We should not forget that, as we 
have already emphasized several times, the energy-momentum tensor for a PhLO must satisfy 
the additional local isotropy (null) condition r^j,(F, *F)r^''(F, *F) = 0. 

So, we have to construct appropriate quantities and relations having direct physical sense in 
terms of the introduced and considered two projections V and V. The above well established in 
electrodynamics relations say that we need two 2-forms to begin with. 

Recall that our coordinate 1-forms dx nd dy have the following vertical and horizontal pro- 
jections: 

V* (dx) = dx — eudz — ud^, H* (dx) = eudz + ud^ , 
V*{dy) = dy — epdz — pd£^, H*{dy) = epdz + pd^. 
We form now the 2-forms V*{dx) A H*{dx) and V*{dy) A H*{dy): 

V*{dx) A H*{dx) = eudx Adz + udx Ad£,, 

V*{dy) A H*{dy) = e p dy A dz + p dx A d^. 
Summing up these last two relations and denoting the sum by F we obtain 

F = eudx Adz + udx Ad£,+ epdy Adz +pdy Ad^. (3.31) 
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Doing the same steps with V* and H* we obtain 

F = —p dx A dz — £ p dx A + u dy A dz + £u dy A d^. (3.32) 

Noting that our definition of the Hodge star requires (*-F)^jy = — ^ £ni/ '^''F^p, it is now easy to 
verify that F = *F. Moreover, introducing the notations 

A = udx+pdy, A* = —epdx + eudy, ( = edz + d^, 

we can represent F and F in the form 

f = aac, f = *f = a* AC. 

From these last relations we see that F and *F are isotropic: F A F = 0, F A *F = 0, i.e. the 
field {F,*F) has zero invariants: F/^i^F^^^ = Fp,v{*FY^ = 0. The following relations are now 
easy to verify: 

V*{F) = H*{F) = V*{*F) = H*{*F) = V*{F) = H*{F) = V*{*F) = H*{*F) = 0, (3.33) 

i.e. F and *F have zero vertical and horizontal projections with respect to V and V. Since, 
obviously, C is horizontal with respect to V and V it is interesting to note that A is vertical 
with respect to V and A* is vertical with respect to V: V*{A) = A, V{A*) = A*. In fact, for 
example, 

V* (A) =V*{udx + p dy) = uV* (dx) + pV* (dy) = 

u[dx + pdz + ep d^] + p[dy — udz — eu d^\ = udx + p dy. 

We are going to establish now that there is real energy-momentum exchange between the 
F-component and the *F-component of the field. To come to this we compute the quantities 
i{Zi)F, i{Z2) * F, i{Zi) * F, i{Z2)F. We obtain: 

i{Z,)F = i{Z2) *F = {A, Z,)C = {A*, Z2)C = \ [{u^ + p\ -e{u^+ p"),] ( = 

= ^F''''{dF),p^dx^ = ^{*Fr{d*F),p^dx^ = iv.T;(F,*F), (3.34) 

i(Zi) * F = -i{Z2)F = {A\ Zi)C = -(A Z2)C = [u{Pi - epz) - p{h -^^^)\ C = 

= -^F"''(d * F),p^dxi' = ]^{*FYP{dF),p^dxi'. (3.35) 

If our field is free then Vi^ T^{F, F) = 0. Moreover, in view of the divergence of the stress-energy- 
momentum tensor given above, these last relations show that some real energy-momentum 
exchange between F and *F takes place: the magnitude of the energy-momentum, transferred 
from F to *F and given by i{Zi) * F = ^(*FYf (dF)ap^dx'^ , is equal to that, transferred from 
*F to F, which is given by —i(Z2)F = — ^^'^''(d * F)fjpf^dx^^. On the other hand, as it is well 
known, the *-invariance of the stress-energy-momentum tensor in case of zero invariants leads 
to F^fjF'^'^ = {*F)ncr{*Fy^ , so, F and *F carry equal and conserved quantities of stress-energy- 
momentum. 

We interpret physically this as follows. The electromagnetic PhLO exist through a special 
internal dynamical equilibrium between the two subsystems of the field, represented by V and 

V, namely, both subsystems carry the same stress-energy-momentum and the mutual energy- 
momentum exchange between them is always in equal quantities. This individualization does 
NOT mean that any of the two subsystems can exist separately, independently on the other. 
Moreover, NO spatial "part" of PhLO should be considered to represent a real physical object. 
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3.6 Electromagnetic PhLO in terms of electromagnetic strain 



The concept of strain is introduced in studying elastic materials subject to external forces of 
different nature: mechanical, electromagnetic, etc. In nonrelativistic continuum physics the 
local representatives of the external forces in this context are usually characterized in terms 
stresses. Since the force means energy-momentum transfer leading to corresponding mutual 
energy-momentum change of the interacting objects, then according to the energy-momentum 
conservation law the material must react somehow to the external interference in accordance 
with its structure and reaction abilities. The classical strain describes mainly the abilities of 
the material to bear force-action from outside through deformation, i.e. through changing its 
shape, or, configuration. The term elastic now means that any two allowed configurations can be 
deformed to each other without appearence of holes and breakings, in particular, if the material 
considered has deformed from configuration Ci to configuration C2 it is able to return smoothly 
to its configuration Ci. 

The general geometrical description [27] starts with the assumption that an elastic material 
is a continuum B C which can smoothly deform inside the space R^, so, it can be endowed 
with differentiable structure, i.e. having an elastic material is formally equivalent to have a 
smooth real 3-dimensional submanifold B C M^. The deformations are considered as smooth 
maps (mostly embeddings) : B ^ M^. The spaces B and are endowed with riemannian 
metrics G and g respectively (and corresponding riemannian co-metrics G""^ and g~^), and 
induced isomorphisms G and g between the corresponding tangent and cotangent spaces . This 
allows to define linear map inside every tangent space of B in the following way: a tangent 
vector V G T^M, re G B, is sent through the differential dip of ip to {dp)x{V) G Ti^^j-^R^, then 
by means of the isomorphism g we determine the corresponding 1-form (i.e. we "lower the 
index"), this 1-form is sent to the dual space r*B of by means of the dual linear map 
{dp))* : T*^^^M.^ — > T*M, and finally, we determine the corresponding tangent vector by means of 

the isomorphism G~^ (i.e. we "raise the index" correspondingly). The so obtained linear map 

:= [G-i o (dip)* o ~g o (dip)] ^ : ^ T^B 

(which is denoted in [27] by (F'^F)^;), extended to the whole B, is called Caushy-Green defor- 
mation tensor field. Now, the combination 

■.= ^[{GoC-G)]^ = ^[{dipro~go{dp)-G]^:TxMxTxM^R 

is called Lagrangian strain tensor field. Note that if we denote by p*g the induced on B metric 
from the metric g (usually euclidean) on then E = ^{p>*g — G). 

We could look at the problem also as follows. The mathematical counterparts of the al- 
lowed (including reversible) deformations are the diffeomorphisms (p oi a riemannian manifold 
{M,g), and every p{M) represents a possible configuration of the material considered. But 
some diffeomorphisms do not lead to deformation (i.e. to shape changes), so, a criterion must 
be introduced to separate those diffeomorphisms which should be considered as essential. For 
such a criterion is chosen the distance change: if the distance between any two fixed points does 
not change during the action of the external force field, then we say that there is no deformation. 
Now, every essential diffeomorphism ip must transform the metric g to some new metric ip*g, 
such that g 7^ p*g. The naturally arising tensor field e = {p*g — 5) 7^ appears as a measure of 
the physical abilities of the material to withstand external force actions. 

Since the external force is assumed to act locally and the material considered gets the cor- 
responding to the external force field final configuration in a smooth way, i.e. passing smoothly 
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through a family of allowed configurations, we may introduce a localization of the above scheme, 
such that the isometry doffeomorphisms to be eliminated. This is done by means of introducing 
1-parametcr group ipt,t G [a, fe] C M of local diffeomorphisms, so, (pa{M) and (pb{M) denote 
correspondingly the initial and final configurations. Now ipt generates a family of metrics (/?| g, 
and a corresponding family of tensors e^. According to the local analysis every local 1-parameter 
group of diffeomorphisms is generated by a vector field on M. Let the vector field X generate 
iff. Then the quantity 

1 ^ 1 ipt Q — g 



2 



ihm^ 
2 t->o t 



i.e. one half of the Lie derivative of g along X, is called (infinitesimal) strain tensor, or defor- 
mation tensor. 

Remark. Further in the paper we shall work with Lx g, i-e. the factor 1/2 will be omitted. 

In our further study wc shall call Lx g, where g = ri is the Minkowski (pscudo)mctric, 
just strain tensor. Clearly, the term "material" is not appropriate for PhLO because no static 
situations are admissible, our objects of interest are of entirely dynamical nature, so 
the corresponding relativistic strain tensors must take care of this. 

According to the preliminary considerations important vector fields in our approach to de- 
scribe electromagnetic PhLO arc ^, A, A*, so, we consider the corresponding three electromag- 
netic strain tensors: L^r]; L^-q; L^, rj. 

Proposition 4. The following relations hold: 



2ux 

Uy+Px 
Uz 



Uy+Px 

2py 

Pz 



Uz 
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{La* ri)^,y = D 





-e{py + Ux) 


-£Pz 


-ep^ 


^{Py + Ux) 


2euy 


euz 




-£Pz 





















Proof. Immediately verified. 

We give now some important from our viewpoint relations. 

Z)(C,C) = D*(C,C) = o, 



D{0 ^ DiO^dx^ 



{u^ - £Uz)dx + {p^ - epz)dy, 



d 
dx>^ 

D 



d 

1 



D^,A*''C 



-{u^ - euz)^ - {p^ - epz)^ = -[A,C], 



^,A'^C = --^{u'+p\-e{u'+p% 



We also have: 



l(p^ - epz) - p{u^ - £u, 

- (Pi - epz)dx + (u^ - euz)dy 
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Clearly, D{Q and D*{Q are linearly independent in general: 

D{Q A D*{() = e {u^- eu^f + {p(_ - sp^f dx Ady = e(j)^{tjj^ - etp^f dx Ady^ 0. 

Recall now that every 2-form F defines a linear map F from 1-forms to 3-forms through the 
exterior product: F{a) := aAF, where a G A^(M). Moreover, the Hodge *-operator, composed 
now with F, gets F{a) back to *F{a) G A^{M). In the previous section we introduced two 
2-forms G = AaC and G* = A* A( and noticed that G* = *G. We readily obtain now 

D{0 ag = d*{Oag* = d{0 aaac = d*{Oaa*ac = 



i{p^ - epz) -p{u^ - euz) 



dx A dy A dz 



t{p^ - epz) - - euz) 



dx Ady Ad^ = 



-((P' L^tp {e dx A dy Adz + dx Ady A d^) 



^ (/y \c / \ uy / \ -r uu. / \ uy / \ uc^ ) — — R {e dx A dy A dz + dx A dy A d^) , 
D{0 AG* = -D*{0 A G = D{0 AA*AC = -D*iO A ^ A C = 

\u'^ + p^)^ - + p'^)z {dx A dy A dz + e dx A dy A dC)- 



Thus, recalling relations (3. 16)- (3. 20), we get 



D{Q A A A C = * D*{C,) A A* A C = -eRC = -i{G*)dG = i{G)dG*, 



D{C) A A* AC 



— * 



D*{C)AAAC = ^ ^ .(^^^^ ^ i{G*)dG*. 



(3.36) 



(3.37) 



The above relations show various dynamical aspects of the energy-momentum redistribu- 
tion during evolution of our PhLO. In particular, equations (3.36-3.37) clearly show that it is 
possible the translational and rotational components of the energy-momentum redistribution 
to be represented in form depending on the (^-directed strains D{Q and D*{C). So, the local 
translational changes of the energy-momentum carried by the two vector components G and 
G* of our PhLO are given by the two 1-forms *[D{() A ^* A C] and *[D*{() A A A C]) and 
the local rotational ones - by the 1-forms *[D{C) A ^ A C] and *[-D*(C) A ^* A C] • In fact, 
the form * D(() A A A (] determines the strain that "leaves" the 2-plane defined by (A, () and 
the form * D*{Q A A* A(] determines the strain that "leaves" the 2-plane defined by (^*,C)- 
Since the PhLO is free, i.e. no energy-momentum is lost or gained from outside, this means 
that the two (null-field) components G and G* exchange locally equal energy-momentum quan- 
tities: * 



D{C) A A AC 



D*{C) A A* AC 



Now, the local energy-momentum conservation 
law [G^o-G''^'^ -|- {G*)^a{G*Y'^~\ = requires L^(j)'^ = 0, and the corresponding strain-fluxes 
become zero: * [D*{C) A A A C] = 0, * [D{C) A A* A C] =0. 

It seems important to note that, only dynamical relation between the local energy-momentum 
change and strain fluxes exists, so NO analog of the assumed in elasticity theory generalized 
Hooke law, (i.e. linear relation between the stress tensor and the strain tensor) seems to exist. 
This clearly goes along with the fully dynamical nature of PhLO, i.e. linear relations exist 
between the divergence terms of our stress tensor ^[ — G^^G^'^ — {G*)fjicr{G*Y'^^ and the C,- 
directed strain fluxes as given by equations (3.36)- (3.37). 
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Chapter 4 

Equations of motion for PhLO. 
Solutions 

In this chapter we show that appropriate solutions for PhLO can be obtained by solving linear 
equations. 

4.1 The approach based on the notion for PhLO 

Every system of equations describing the time-evolution of some physical system should be con- 
sistent with the very system in the sense that all identification characteristics of the system 
described must not change. In the case of electromagnetic PhLO we assume the couple {F, F) 
to represent the field, and in accordance with our notion for PhLO one of the identification 
characteristics is straight-line translational propagation of the energy-density with constant ve- 
locity "c", therefore, with every PhLO we may associate appropriate direction, i.e. a geodesic 
null vector field = on the Minkowski space-time. On the other hand, the complex of 
field functions {Ffj_i^ , Ffj_i^) admits both translational and rotational components of propagation. 
We choose further ^ = —£-§^ + which means that we have chosen the coordinate system in 
such a way that the translational propagation is parallel to the plane iz,^,). For another such 
parameter we assume that the finite longitudinal extension of any PhLO is fixed and is given by 
an appropriate positive number A. In accordance with the "compatible translational-rotational 
dynamical structure" of PhLO we shall assume that no translation is possible without rotation, 
and no rotation is possible without translation, and in view of the constancy of the translational 
component of propagation we shall assume that the rotational component of propagation is 
periodic, i.e. it is characterized by a constant frequency. The natural period T suggested is ob- 
viously T = ^. An obvious candidate for "rotational operator" is the linear map J transforming 
F to F, which map coincides with the reduced to 2-forms Hodge-*. Geometrically, * rotates the 
2-frame {A, A*) to |, so if such a rotation is associated with a translational advancement of to, 
then a full rotation should correspond to translational advancement of ilo = A. The simplest 
and most natural translational change of the field (F, F) along ^ should be given by the Lie 
derivative of the field along ^. Hence, the simplest and most natural equations should read 

KloL^{F)=sF, (4.1) 

where F and F are given in the preceding section, k = ±1 is responsible for left/right orientation 
of the rotational component of propagation, and lo = const. Vice versa, since J o J = —id and 
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J ^ = — J the above equation is equivalent to 

kIo L^{F) = -eF. 

It is easy to show that these equations are equivalent to 

kIo L^{V - Vo) = e{V - Vo), (4.2) 

where V is given by (40) and Vg = dx + dy ^ in our coordinates is the identity map in 
Im{V) = Im{V). Another equivalent form is given by 

kIoZi = a*, or KI0Z2 = — A, 

where A* and A are ry-corresponding vector fields to the 1-forms A* and A. 

4.2 The Lagrangian Approach 

Appropriate lagrangian for the above equations (/o=const.) is 

= i (^'"""^ - ^-^) - i + ^"'') 

F and F are considered as independent, and the relations FAF = FAF = FAF = lead to 

^g dFap ^ -^dFa^^^p ^ ^ 

^ dx" ^ dx" 

The corresponding Lagrange equations read 



^«^or ^ + Fap = 0, enloC^-^ - F^, = 0, (4.4) 



so, on the solutions the lagrangian becomes zero: 'L,{solutions) = 0. The stress-energy-momcntum 
tensor, in view of the null character of F and F, is given by (3.29), where *F has to be replaced 
by F. It deserves noting that the above null conditions lead to F^„F'^'^ = F/j^F^'^ and to 
F^(jF'^'^ = 0. Hence, the two subsystems represented by F and F carry the same strcss-cncrgy- 
momentum, therefore, F F cncrgy-momentTim exchange is possible only in equal quantities. 
In our coordinates the above equations reduce to 

kIo{u^- euz) = -p, kIo{p^ - e pz) = u, 

it is seen that the constant Ip satisfies the above given relation (3.28). Prom these last equations 
we readily obtain the relations 

(n^ +/)e -eiu^ +p^)z = 0, u{p^ - epz)-p{u^-euz) = ^{u^ +p^), 

which represent our equations in energy-momentum terms. Now, the substitution u = <t cos ijj, p = 
^sintp, leads to the relations 

Lf$ = 0, LcV' = f- 

In terms of nonlinear connections recalling that = —^tr{V o H*) and computing ^tr{V o 
L^H*) = £[u{p^ — epz) — p{u^ — £Uz)] = ^'^eL^ijj the last two relations can be equivalently 
written as 

[tr{V o H*)] = 0, tr{V o L^H*) = -^tr{V o H*). 
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4.3 Equations of motion in terms of translational-rotational 
compatability 

In order to look at the translational-rotational compatability as a generating tool for writing 
equations of motion we recall first the concept of local symmetry of a distribution: a vector field 
y is a local (or infinitesimal) symmetry of a p-dimensional distribution A defined by the vector 
fields {Yi,...,Yp) if every Lie bracket is in A: [Fj,!^] € A. Clearly, if A is completely 

integrable, then every is a symmetry of A, and the flows of these vector fields move the points 
of each completely integral manifold of A inside this integral manifold, that's why they are called 
sometimes internal symmetries. If Y is outside A then it is called shuffling symmetry [28], and 
in such a case the flow of Y transforms a given completely integral manifold to another one, i.e. 
the flow of Y " shuffles" the lists of the corresponding foliation. We are going to show that our 
vector field ( = —e-^ + is & shuffling symmetry for the distribution Aq deflned by the vector 

fields {A, A*). In fact, Ag coincides with our vertical distribution generated by (^, so it is 
completely integrable and its integral manifolds coincide with the (x, ?/)-planes. From physical 
point of view this should be expected in view of the intrinsically required stability of our PhLO 
under translational propagation along null straight lines: this propagation just transforms the 
2-plane {x,y) passing through the point to a parallel to it 2-plane passing through the 

point (2:25^2), and these two points lay on the same trajectory of our field ^. 
The corresponding Lie brackets are 



d 



d - - d d 



We see that [A, (] and [^*, C] are generated by but C is outside Ag, so our field C is a 

shuffling local symmetry of Aq. 

We notice now that at each point we have two frames: {A, A* , dz, d^) and ([A, Q, [A*,(],dz,d^). 
Since physically we have internal energy-momentum redistribution during propagation, we could 
interprete the permanent availability of these two intrinsically connected frames as correspond- 
ing mathematical adequate of this physical process. Taking into account that only the first two 
vectors of these two frames change during propagation we write down the corresponding linear 
transformation as follows: 



{[A,C],[A*,C]) = {A,A* 



Solving this system with respect to (a, /3, 7, 5) we obtain 



a (3 


1 


7 6 


~ 02 



-£R 



eR 



'2~52~ 



■1 



where R = n (p^ — epz) — p{u^ — £Uz)- If the translational propagation is governed by the 
conservation law L^^^ = 0, then we obtain that the rotational component of propagation is 
governed by the matrix eL^'tpJ, where J denotes the canonical complex structure in M?, and 
since L^ijj = u{p^ — epz) — p{u^ — suz) 7^ we conclude that the rotational component 
of propagation would be available if and only if R 7^ 0. We may also say that a compatible 
translational-rotational dynamical structure is available if the amplitude <I>^ = + p^ is a 
running wave along C, and the phase ip = arctg^ is NOT a running wave along C, : L^^i}) 7^ 0. 
Physically this means that the rotational component of propagation is entirely determined by 
the available internal energy-momentum exchange: i{F)dF = —i{F)dF. 



50 



Now, if we have to guarantee the conservative and constant character of the rotational aspect 
of the PhLO nature, we can assume L^ip = const = kI~^,k = ±1. Thus, the frame rotation 
{A,A*,d„di:) {[A,C],[A*,C],d„d^), i.e. [A,C] = -eA*L-^iP and [^*,C] = eAL-^^p, gives the 
following equations for the two functions (n,p): 

K K 

u^- euz = - — p, p^ - epz = Y u . 

The quantity R = n (p^ — epz) —p{u^ — e Uz) = ^'^L^tp = kI^^^^ suggests to find an integral 
characteristic of the PhLO rotational nature. In fact, the two co-distributions {A, ^) and {A*,Q 
define the two (equal in our case) Probenius 4-forms dA A A A C = dA* A ^* A C- Each of these 
two 4-forms is equal to eRwo- Now, multiplying by lo/c any of them we obtain: 

-dAAAAC = -dA*AA*AC = -eBxJo = sk—Uo . (4.5) 
c c c c 

Integrating over the 4-volume x (A = AIq) (and having in view the spatially finite nature of 

PhLO) we obtain the finite quantity TC = £kET, where E is the integral energy of the PhLO, 

r = ^, which clearly is the analog of the Planck formula E = hu, i.e. h = ET. The combination 

£K means that the two orientations of the rotation, defined by k = ±1, may be observed in each 

of the two spatial directions of translational propagation of the PhLO along the z-axis: from 

—DO to +00, or from +00 to —00. 

Finally, recalling relations (3.15), we can easily see that in case of L^ip = kI~^ and = 

the 3-form SF A F = eR * ( = * ( becomes closed: d{dF A F) = 0, which also gives an 

integral conservation law. In fact, the 3-integral of the reduced on 3-form -^{SF A F) gives 
£kET, where E is the integral energy, so, the Planck formula holds. 

4.4 Photon-like Solutions 

We consider the equations obtained in terms of the two functions $ = \/v? -\-p^ and V' = arctg^. 

The equation for $ in our coordinates is $g — e^z = 0, therefore, $ = ^{x, y,(,+ £z), where $ is 
allowed to be spatially finite, as assumed further, or spatially localized function. The equation 
for 'tp is tp^ — £ipz = Two families of solutions for tp, depending on an arbitrary function </? 
can be given by 

tpi = -^z + (fix, y,^ + £z), and 1^2 = + (pix,y,^ + ez). 

Since is a spatially finite function representing the energy density we see that the translational 
propagation of our PhLO is represented by a spatially finite running wave along the z-coordinate. 
Let's assume that the phase is given by V"! and, for simplicity, = 0. The form of this solution 
suggests to choose the initial condition ut=o{x,y,£z),pt=o{x,y,£z) in the following way. Let for 
z = the initial condition be located on a disk D = D{x,y;a,b;ro) of small radius To, the 
center of the disk to have coordinates (a, 6), and the value of ^t=Q(x,y,0) = «|^o + Pt=o 
be proi)ortional to some appropriate for the case bump function / > on D of the distance 
(x — a,y^ + {y — between the origin of the coordinate system and the point (a;,y, 0), such 
that it is centered at the point (a, 6), so, f{x,y) = f{\/{x — a)^ + (y — fc)^ ), D is defined by 
D = {(x, y)\\/ (x — a)2 + {y — by < To}, and f{x,y) is zero outside D. Let also the dependence 
of $t=o on z be given by be the corresponding bump function 9{z; A) > of an interval {z, z + X) 
of length A = 4lo on the z-axis. If 7 > is the proportionality coefficient we obtain 

= 7 /(x, y; a, b) 9{ct + £z; A) cos('i/'i), 
p = 7 f{x, y; a, b) e{ct + £z; A) sin(V'i). 
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We see that because of the available sine and cosine factors in the solution, the initial condition 
for the solution will occupy a Sd-spatial region of shape that is close to a helical cylinder of 
height A, having internal radius of and wrapped up around the z-axis. Also, its center will 
always be \/a^~-j-~P^-distant from the z-axis. Hence, the solution will propagate translationally 
along the coordinate z with the velocity c, and, rotationally, inside the corresponding infinitely 
long helical cylinder because of the z-dependence of the available periodical multiples. 

On the two figures below are given two theoretical examples with k = —1 and k = 1 
respectively, amplitude function $ located inside a one-step helical cylinder P with height of A, 
and phase ip = —en'^^ = —sk^. The solutions with e = — 1 will propagate left-to-right along 
the coordinate z. 





Figure 4.2: Theoretical example with k = 1. The translational propagation is directed left-to- 
right. 



The curvature K and the torsion r of the screwline inside D through the point (x, y,0) £ D will 
be 

(7/0)2 + 62' ^ (^/^)2 + 52 ' 

where b = X/2tt = IIq/tt. The rotational frequency will he v = c/X = c/AIq, so we can 
introduce period T = l/v and elementary action h = E.T, where E is the (obviously finite) 
integral energy of the solution defined as 3d-integral of the energy density $2 = (7/^)^. 
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Chapter 5 



Retrospect 



In trying to understand our observational knowledge of the real world we must be able to separate 
the important structural and behavioral properties of the real objects from those, the changes 
of which during time-evolution do not lead to annihilation of the objects under consideration. 
One of the basic in our view lessons that we more or less have been taught is that the physical 
objects are spatially finite entities, and that for their detection and further study some energy- 
momentum exchange is necessarily required. So, every physical object necessarily carries energy- 
momentum and every interaction between two physical objects has such an energy-momentum 
exchange aspect. The second lesson concerning any interaction is that, beyond its universality, 
energy-momentum is a conserved quantity, so NO loss of it is allowed: it may only pass from 
one object to another. This means that every annihilation process causes creation process(es), 
and the full energy-momentum that has been carried by the annihilated objects, is carried 
away by the created ones. Energy-momentum always needs carriers, as well as every physical 
object always carries energy-momentum. Hence, the energy-momentum exchange abilities of any 
physical object provide protection against external influence on one side, and reveal its intrinsic 
nature, on the other side. Therefore, our knowledge about the entire complex of properties 
of a physical object relies on getting information about its abilities in this respect and finding 
corresponding quantities describing quantitatively these abilities. 

The spatially finite nature of a physical object implies spatial structure and finite quan- 
tity of energy- momentum needed for its creation, so NO structureless and infinite objects may 
exist. The approximations for "point object" and "spatially infinite field", although useful in 
some respects, seem theoretically inadequate and should not be considered as basic ones. More 
reliable appears to be the approximation "finite continuous object", which we tried to follow 
throughout our exposition. This last approximation suggests that the usually nonlinear local 
energy-momentum conservation relations should be the basic suggesting tool for finding ap- 
propriate mathematical models of local nature of such objects. The natural physical sense of 
the corresponding model equations is not necessarily supposed to be local energy-momentum 
exchange, but must be consistent with it, and closely connected to it. 

Another useful observation is that physical objects are many-aspect entities, they have com- 
plicated structure and their very existence is connected with internal energy-momentum ex- 
change among the various structural components. So, the mathematical model objects should 
be many-component ones, and with appropriate mathematical structure. Of basic help in find- 
ing appropriate mathematical objects is having knowledge of the internal symmetry properties 
of the physical object under consideration. This "step by step" process of getting and accu- 
mulating important information about the physical properties of natural objects reflects in the 
"step by step" process of refining the corresponding mathematical models. 
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The greatest discovery at the very beginning of the last century was that the notion of 
electromagnetic field as suggested by Maxwell equations is inadequate: the time dependent elec- 
tromagnetic field is not an infinite smooth perturbation of the aether, on the contrary, it consists 
of many individual time-stable objects, called later photons, which are created/destroyed mainly 
during intra-atomic energy-transition processes. Photons are finite objects, they carry energy- 
momentum and after they have been radiated outside their atom-creator, they propagate as a 
whole translationally by the speed of light. Moreover, their propagation is not just translational, 
it includes rotational component, which is of intrinsic and periodical nature. The corresponding 
intrinsic action for one period T is h = ET, where E is the full energy of the photon, and all 
photons carry the same elementary intrinsic action h. During the entire 20th century physi- 
cists have tried to understand the dynamical structure/nature of photons from various points 
of view, and this process is still going on today. The conviction that a new point of view on 
the dynamical nature of the field equations is needed is shortly summarized by Ziolkowsi [29]: 
"Finite-energy, diffraction- free beams for the linear free-space wave equation are imposible". 

In order to come to a new look at the situation in nonrelativistic terms we made use of the 
Newton approach: the identifying features of the object considered must be kept unchanged 
during evolution and the admissible changes most naturally should be expressed by means of 
specializing the energy-momentum exchange abilities of the object considered, as well as by 
paying due respect to the available translation-rotation interrelation. We showed that this 
approach works well in the nonrelativistic (E, B) formalism and concluded that the relativistic 
{F, *F) structure is much more adequate to the PhLO dynamical structure than the (E, B) one. 

The basic theoretical idea in the relativistic formalism was to make use of the Frobenius inte- 
grability/nonintegrability theorems as an appropriate mathematical machinary: the integrability 
of a distribution we connected with the time-stability of the basic identification properties of 
the object considered, and the nonintegrability of the available subdistributions was interpreted 
physically as internal interaction among the subsystems, where the corresponding curvatures 
appeared as natural mathematical tools for generating appropriate mathematical images of the 
local energy-momentum exchange fluxes between any two subsystems. 

We introduced a notion of PhLO as a spatially finite physical object with a compatible 
translational-rotational dynamical structure and propagating translationally with the frame in- 
dependent velocity of light "c". We showed that Frobenius integrability theory possesses all 
necessary features to meet the physical aspects of this notion. From physical viewpoint, two 
dynamically interacting subsytems of a PhLO can be individualized, these subsystems carry the 
same strcss-cncrgy-momcntum, and they exchange energy-momentum locally always in equal 
quantities, so they exist in a dynamical equilibrium,. The mathematical realization of the two 
subsystems of a PhLO was made in two ways: through a direct choice of two nonintegrable sub- 
distributions, and by means of a couple of two nonlinear connections V and V with a common 
image space. Their inter-communication is carried out and guaranteed by the nonzero curvature 
forms and 17* in the first case, and by the nonzero curvature forms TZ and TZ in the second 
case. The values of these curvature forms define two 1-dimensional space-like subspaces, so, the 
corresponding two exterior products with the null direction of translational propagation give the 
mathematical images F and F of the two interacting subsystems. This approach allows to get 
some information concerning the dynamical nature of the PhLO structure not only algebraically, 
but also infinitesimally, i.e. through the curvature forms. 

While the energy density of a PhLO propagates only translationally along straight 
isotropic lines, the available interaction of the two subsytems of a PhLO demonstrates itself 
through a rotational component of the entire propagational behaviour and is available only if 
the curvature forms are not zero. The mutual energy-momentum exchanges are given by the 
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inner products of the curvature images with F and *F. The dynamical equihbrium between F 
and *F is given by i{*F){dF) = —i{F){d * F). 

Besides the spatially finite nature of PhLO that is allowed by our model and illustrated 
with the invariant parameter Ig, two basic identifying properties of PhLO were substantially 
used: straight-line translational propagation with constant speed, and constant character of the 
rotational component of propagation. The physical characteristics of a PhLO are represented by 
an analog of the Maxwell-Minkowski stress-energy-momentum tensor. An interesting moment 
is that F and F have zero horizonal and vertical components with respect to the two nonlinear 
connections. 

It was very interesting to find that some of the basic characteristics of PhLO could be given 
in terms of the two strain-tensors, i.e. through the Lie derivatives of the Minkowski pseudo- 
metric with respect to spatial direction generators of the two nonintegrable subdistributions, so, 
each of the two nonintegrable subdistributions has its own strain tensor. This gives entirely new 
viewpoint on PhLO, namely, the PhLO energy-momentum propagates through deformations! 
The values of each strain tensor on the generators of its subdistribution gives the translational 
change of the energy density, and its value on the generators of the other subdistribution gives 
the Probenius curvature. The corresponding dynamical aspects are given by expressions (3.36)- 
(3.37). 

It seems important to note that the curvature forms are not zero only if the component- 
functions of the vector fields defining the distributions (or the component-functions of the as- 
sociated nonlinear connections) are NOT running waves along the translational propagation, 
e.g. the squares of the curvature forms are equal to |dFp = |d * F\'^ = \5F\'^ = \5 * F\'^ = 
{u^ — £Uz)^ -\- ipi,— ^Pz)^ 7^ 0. Also, the dually invariant longitudinal size parameter acquires 
sense only for finite nonzero curvatures. 

The equations of motion can be viewed from different viewpoints: as compatability condi- 
tions between the rotational and translational components of propagation, as Lagrange equa- 
tions for an action principle, as the nonlinear part of the solutions of the vacuum equations 
of EED, and also as naturally defined transformation of 2-dimensional frames. In all these 
aspects of the equations of motion the curvature forms play essential role through controlling 
the inter-communication between F and F. Moreover, the Frobenius curvature turns out to be 
proportional to the energy density, which recalls the main idea of General Relativity from one 
side, and allows an analog of the famous Planck formula E = hf io be introduced, from the 
other side. 

The solutions considered illustrate quite well the positive aspects of our approach. It is 
interesting to note that the phase terms of these solutions depend substantially only on spa- 
tial variables, so, the spatial structure of the solutions considered participates directly in the 
rotational component of the PhLO dynamical structure. 

Our basic conclusion reads: PhLO are complex objects with dynamical structure of special 
kind, so any mathematical model of PhLO shall need corresponding mathematical structure. 

According to the results given in this study the basic adequate mathematical structure in case 
of electromagnetic PhLO is the 3-dimensional distribution/codistribution {A, A* ,C,) / {A, A* ^C,) 
on Minkowski space-time together with the corresponding interconnections represented by the 
integr ability /nonintegr ability properties of its subdistributions. 

This study was partially supported by Contract ^1515/2005 with the Bulgarian National 
Fund " Science Research" . 
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